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Preface 


I! 


The  Electronic  Warfare  Division  (WR)  of  the  Air  Force 
Avionics  Laboratory  expressed  a need  for  a model  to  predict 
the  atmospheric  effects  on  a slightly  divergent  laser  beam 
which  is  reflected  from  a small  comer*cube  and  received  bach 
at  the  receiver.  Such  an  experiment  was  being  set  up,  with 
beam  parameters  as  described  in  this  report.  I undertook  the 
task  of  developing  a model  adequate  for  describing  the  exper- 
iment and  hopefully  other  cases  of  interest. 

The  study  of  atmospheric  phenomena  was  a completely  new 
one  for  me,  and  you  will  find  no  original  contributions  here. 
This  study,  and  this  model,  are  basically  a compilation  of 
other  peoples  results,  as  applied  to  a particular  problem. 

A list  of  people  to  whom  1 am  indebted  for  help  and  guid- 
ance would  have  to  start  with  my  entire  bibliography. 

On  a more  personal  level,  however,  my  thanks  go  toi 
Major  Bruce  Pierce,  my  thesis  advisor,  for  his  suppoirt  and 
guidance,  and  Major  Ron  Grotbeck,  of  the  Avionics  Lab,  for 
his  help  with  the  description  of  the  experimental  system  and 
with  finances, 

1 am  also  Indebted  to  Professor  Donn  Shankland  for 
attempting  to  guide  me  through  Tatarskl ' s labsrrlnthean 
work,  and  to  S.  A.  (Tony)  Clough,  and  many  others  at  the  Air 
Force  Geophysics  Laboratory,  for  trying  to  show  me.  In  one 
week,  how  to  handle  the  atmospheric  molecular  absorption 
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aodels  which  they  have  taken  years  to  develop. 

On  a still  more  personal  level*  I want  to  thank  my 

wife*  Marilyn,  and  our  children,  Kim  and  Scott,  for  their 

tolerance  and  understanding  over  the  past  few  months. 

William  H.  Macinnis 
Major  USAF 
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Abstract 

A model  is  described  which  calculates  the  irradiance 
probability  distribution  at  the  receiver  of  a system  con- 
sisting of  a co-located  laser  transmitter  and  receiver, 
with  the  beam  being  reflected  from  a small,  highly  reflective 
target.  The  beams  considered  have  short  pulses  (10  nsec)  at 
low  rates  (30  Hz),  Atmospheric  effects  included  are  molecular 
scattering  and  absorption,  aerosol  extinction  and  turbulence. 
It  is  found  that,  in  normal  daytime  weather  conditions,  it 
is  common  to  have  pulse -to -pulse  irradiance  variations  at 

< 

the  receiver  of  two  orders  of  magnitude. 

! 

i 
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I,  Introduction 


Background 

The  Electronic  Warfare  Division  of  the  Air  Force  Avi- 
onics  Laboratory  is  currently  conducting  experimental 
studies  of  pulsed  laser  beams  to  determine  the  feasibility 
of  using  such  beams  in  future  countermeasures  applications. 
The  basic  series  of  experiments  entails  measuring  the  effects 
of  the  atmosphere  on  pulsed  beams  over  a test  range  at 
Wright -Patterson  Air  Force  Base,  Ohio. 

The  laser  being  used  operates  in  the  visible  and  near- 
infrared  wavelength  region.  It  emits  beams  at  the  Neod3Tniura 
wavelength  near  1.06  microns,  the  frequency  doubled  wave- 
length near  .53  microns^ and  a variable  wavelength  between 
the  other  two.  The  variable  wavelength  beam  comes  from  a 
tuneable  dye  laser,  pumped  by  the  Nd  beam.  Each  beam  is 
slightly  divergent  (several  milliradians).  The  beams  are 
reflected  from  a small  comer  reflector  at  the  other  end  of 
the  range  and  detected  by  a receiver  collocated  with  the 
transmitter.  These  are  low  duty-cycle  pulsed  beams  with 
pulse  durations  of  less  than  ten  nanoseconds  and  pulse  repe- 
tition frequencies  of  about  30  Hz.  Hence,  peak  powers  on 
the  order  of  ten  megawatts  are  developed  with  average  powers 
of  only  about  three  watts. 
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There  has  been  no  model  of  atmospheric  propagation 
which  fully  describes  the  propagation  of  such  beams.  Most 
recent  effort  has  been  directed  toward  Investigating  con- 
tinuous wave  or  high  duty-cycle  pulsed  lasers  for  energy 
transfer  (Including  weapons)  or  communications.  In  partlc- 
ular»  work  has  been  done  on  compensating  for  the  random 
effects  of  atmospheric  turbulence  In  the  optical  communi- 
cations field  and  energy  losses  and  defocuslng  effects  have 
been  widely  considered  for  energy  transfer.  The  AGARD  Con- 
ference proceedings  (Ref  1)  from  1976  descrll}e  the  then  cur- 
rent state  of  Investigations  Into  various  aspects  of  the 
optical  propagation  problem. 

The  various  phenomena  which  have  to  be  taken  Into  2iccount 
In  order  to  describe  optical  beam  propagation  through  the 
atmosphere  arei 

1.  scattering  and  absorption  by  atmospheric  gases 
(called  molecular  scattering  and  absorption) 

2.  scattering  and  absorption  by  suspended  solid  or 
liquid  particles 

(collectively  called  aerosol  extinction) 

3.  scattering  due  to  random  variations  In  the  local 
Index  of  refraction  of  the  atmosi^ere 

(called  tvurbulence) . 

It  should  be  noted  that  turbulence.  In  this  context, 
refers  to  random  Inhomogeneltles  In  the  refractive  Index, 
caused  primarily  by  corresponding  Inhomogeneltles  In  the 
local  temperature  of  the  air  mass.  This  Is  not  the  same  as 
what  Is  normally  called  turbulence)  l.e. , random  fluctuations 


In  the  local  wind  velocity.  Indeed,  such  wind  flow  irregu- 
larities may  well  provide  good  mixing  and  thus  reduce  temper- 
atijre  inhomogeneities. 

Each  of  the  factors  above  has  been  widely  studied.  A 
recent  international  conference,  mentioned  above  (Ref  1), 
was  devoted  to  sharing  recent  results  and  discussing  needed 
efforts  throughout  the  field  of  optical  propagation.  There 
is  active  research  going  on  in  all  of  the  above  areas,  but 
the  degree  of  understanding  and  the  data  base  against  which 
hypotheses  can  be  checked  varies  considerably  among  them. 

Ob lectives 

The  goal  then,  is  to  develop  a model  which  will  predict 
on  a pulse-to-pulse  basis,  the  power  at  the  receiver,  based 
upon  normally  reported  and  forecast  meteorological  variables. 
For  any  application  relying  upon  reflected  power  measurements 
for  tracking  or  aiming,  it  is  important  to  know  what  sort  of 
variation,  pulse-to-pulse,  can  be  expected  solely  on  the 
basis  of  the  atmospheric  effects.  Determination  of  the 
acceptable  level  of  variation  is  not  a part  of  this  effort. 
That  is  clearly  a function  of  the  nature  of  the  reflector, 
the  installation  of  the  laser  device  and  the  tracking  algo- 
rithm involved  in  a real  system.  For  the  present  fixed 
installation  with  a stationary  comer  reflector,  tracking  is 
not  a problem  and  the  aiming  is  considered  to  be  optimum. 

The  model  must  take  into  account  all  of  the  mechanisms 
described  above.  It  must  give  both  the  average  power  received 
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and  the  probability  distribution  of  the  power.  Experimental 
results  are  plotted  as  the  number  of  pulses  in  various  power 
intervals,  after  counting  a large  number  of  pulses.  Hence, 
these  results  should  essentially  have  the  same  form  as  the 
probability  distribution,  if  the  intervals  are  suitably 
chosen. 

Approach 

Each  of  the  above  factors  is  considered  separately.  Fol- 
lowing the  usual  approach  (Ref  23),  it  is  assumed  that  Beer's 
Law  is  appropriate  for  the  description  of  attenuation  due  to 
molecular  and  aerosol  effects.  Hence,  the  irradiance  reaching 
the  plane  at  a distance  x from  the  source  (of  initial  irradi- 

ance  E^)  is 
eo 

Ee  " ^eo  ® 

and 

- o.  + '*■  ^a 

idiere 

is  the  molecular  scattering  coefficient 
is  the  molecular  absorption  coefficient 

and 

is  the  combined  aerosol  extinction  coefficient. 

Each  of  the  three  coefficients  is  calculated  by  use  of 
information  which  is  the  most  current  available.  Molecular 
scattering  is  based  upon  Rayleigh  scattering  theory,  molecular 
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absorption  upon  the  Air  Force  Geophysics  Laboratory  (AFGL) 
Absorption  Line  Compilation  (Ref  21)  and  aerosol  extinction 
upon  the  latest  aerosol  models  incorporated  into  LOWTRAN  3B 


(Refs  30,  31). 
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It  is  then  assumed  that  the  molecular  and  aerosol  effects 
are  not  coupled  to  the  turbulence  effects.  That  is,  it  is 
assumed  that  the  total  atmospheric  effect  is  the  same  as 
though  the  power  were  reduced  by  the  molecular  and  aerosol 
attenuation  first,  and  then  the  beams  were  subjected  to  the 
turbulent  path.  The  various  effects  of  turbulence  upon  the 
beam  are  investigated  and  the  degree  to  which  each  affects 
this  particular  problem  is  discussed.  Those  of  importance 
are  then  calculated  and  applied  to  the  determination  of  what 
radiation  is  incident  upon  the  target. 

The  target  reflectivity  and  divergence  are  then  applied 
and  the  reflected  beam  is  essentially  treated  with  another 
whole  application  of  the  same  procedure. 

The  procedures  and  formulas  used  for  turbulence  calcu- 
lations are  not  nearly  as  clear  cut  as  those  for  the  other 
phenomena.  In  particular,  in  moderately  strong  turbulence, 
there  is  no  general  agreement  even  as  to  the  nattuce  of  the 
irradiance  probability  distribution.  One  of  the  proposed 
distributions  (Ref  9)  is  chosen,  based  upon  agreement  with 
closely  controlled  experiments.  Throughout  the  rest  of  the 
turbulence  calculations,  approaches  which  have  apparently 
given  adequate  results  in  similar  applications  are  used. 
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In  all  the  calculations,  no  attempt  is  made  to  derive 
the  formulas  or  to  develop  new  theories.  The  approach, 
rather,  is  strictly  a pragmatic  one  of  selecting  suitable 
results  from  other  sources  and  applying  them  to  this  par- 
ticular problem. 

The  computer  program  is  meeint  to  be  simple  to  use  and 
self-explanatory  in  the  output,  to  give  a good  idea  of  just 
%fhat  can  be  expected  in  system  performance.  The  program  is 
designed  on  a strictly  modular  basis,  with  all  calculations 
done  in  subroutines.  Hence,  if  it  is  later  determined  that 
more  accuracy  is  required  in  some  area,  or  if  new  approaches 
are  discovered,  only  the  appropriate  subroutines  need  be 
changed.  For  example  (Ref  24) , work  is  currently  under  way 
at  AFGL  on  improved  aerosol  models  incorporating  a dependence 
on  relative  humidity.  This  could  be  incorporated  with  a 
minimum  of  trouble. 
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The  propagation  of  electromagnetic  radiation  through 
the  atmosphere  is  a highly  complex  subject,  some  aspects  of 
which  are  not  well  understood.  There  are  several  different 
mechanisms  involved,  which  vary  in  significance  depending 
upon  both  the  nature  of  the  radiation  and  the  nature  of  the 
propagation  path  being  considered. 

The  radiation  being  considered  in  the  present  case  is 
optical  in  nature,  with  wavelengths  in  the  visible  and  near- 
infrared.  It  is  in  the  form  of  pulsed  laser  beams  which  are 
slightly  diverging.  The  beams  consist  of  pulses  of  duration 
(PD)  on  the  order  of  1-10  nanoseconds  and  with  a repetition 
rate  of  about  30  Hz, 


The  effects  which  are  of  importance  can  be  grouped  into 
three  main  areas t molecular  effects,  aerosol  effects^  and 
refractive  effects. 

Molecular  Effects 

In  the  context  of  atmospheric  physics,  the  term  "molec- 
ular” generally  refers  to  the  gas  molecules  making  up  the 
atmosphere,  as  opposed  to  any  suspended  particulates,  called 
aerosols.  Beams  propagating  through  the  atmosphere  suffer 
losses  due  to  two  different  molecular  mechanisms --scattering 
and  absorption. 
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Molecular  Scatter.  The  scattering  of  electromagnetic 
radiation  particles,  small  compared  to  the  wavelength  of 


the  radiation,  is  generally  called  Rayleigh  scattering.  Such 
scattering  has  been  extensively  studied  (Ref  4t651>652)  and 
results  are  available  concerning  all  aspects  of  the  effect. 
Since  only  the  amount  oif  light  scattered  is  of  importance 
here,  the  interesting  angular  and  polarization  effects  are 
not  considered.  In  the  wavelength  range  of  interest,  the 
attenuation  due  to  molecular  scattering  losses  is  very  small, 
on  the  order  of  10  per  kilometer t i.e.,  over  one  kilometer, 
the  transmission  of  a beam  will  be  exp(>10  ) or  99.9%,  when 

only  molecular  scatter  is  considered.  The  molecular  scatter 
effect  is  included  for  completeness  even  though  the  effect 
is  negligible,  except  tinder  the  very  clearest  conditions 
when  aerosols  are  not  a factor. 

Molecular  Absorption.  Unlike  molecular  scatter,  tdiich 
is  a smooth,  slowly  varying  function  of  wavelength,  molecular 
absorption  shows  very  strong  wavelength  dependence.  This  is 
because  it  is  a qtiantum  effect,  in  which  each  molecular 
species  present  can  absorb  only  radiation  of  energies  corre- 
sponding to  its  own  internal  energy  level  differences. 

The  amount  of  absorption,  then,  clearly  depends  not 
only  upon  wavelength  but  also  upon  the  amount  of  each 
absorbing  species  present  in  the  path.  The  Air  Force  Geo- 
physical Laboratory  has  been  compiling  data  on  both  the 
makeup  of  the  atmosphere  and  the  line  parameters  for  the 
absorption  lines  of  each  species  (Refs  21,  22).  This  latter 
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compilation  has  drawn  from  both  laboratory  measurements  of 
absorption  and  from  quantum  mechanical  calculations.  Chap- 
ter III  covers  in  detail  the  calculation  of  the  attenuation 
from  this  data  base. 

Aerosol  Attenuation 

The  effects  of  aerosol  attenuation  are  those  with  which 
everyone  is  familiar.  Haze,  fog,  smoke,  and  clouds  all  fall 
into  this  category,  hence  its  importance  to  propagation  is 
obvious.  Methods  of  accounting  for  the  effect,  however,  are 
much  less  obvious.  In  a paper  published  in  1908  (Ref  4i634), 
Mie  obtained  from  basic  electromagnetic  theory  a rigorous 
solution  for  diffraction  of  a monochromatic  plane  wave  by  a 
homogeneous  sphere  of  any  size  and  composition,  suspended  in 
a homogeneous  medium.  Although  most  atmospheric  aerosol 
particles  are  not  spherical,  the  calculations  based  on  this 
solution  yield  good  information  concerning  the  aerosol  prob- 
lem, In  fact,  it  has  been  shown  (Ref  33il9)that  extinction  by 
non- sjAierical  aerosol  particles  is  less  than  that  for  spher- 
ical particles,  so  that  use  of  Mie's  theory  yields  a con- 
servative estimate  of  transmittance. 

The  basic  Mie  formulation  is  for  a single  particle,  but 
it  can  be  extended  to  consideration  of  a distribution  of 
particles,  so  long  as  they  are  randomly  distributed  and 
separated  by  distances  large  compared  to  a wavelength 
(Ref  4 i634).  Hence,  it  is  possible  to  calculate  quite 
rigorously  both  the  scattered  and  absorbed  radiation,  given 
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the  aerosol  distributions.  The  calculations  themselves  are 
far  from  trivial,  since  even  the  relatively  simple  case  of  a 
single  particle  results  in  a lengthy  solution  in  terms  of 
infinite  series,  which  are  not  always  rapidly  convergent, 

A more  basic  problem  in  the  aerosol  calculations  is  the 
need  to  know  the  distribution  of  particles,  A rigorous 
solution  requires  knowledge  of  the  size,  refractive  index 
and  number  density  of  particles  as  a function  of  position. 

The  difficulties  found  in  determining  these  parameters  and 
the  models  used  in  this  study  are  found  in  Chapter  IV, 

Refractive  Effects 

The  bulk  properties  of  the  atmosphere,  as  described  by 
the  refractive  index,  also  have  an  effect  on  the  propagation. 
There  are  both  smoothly  varying  and  random  effects  present. 

Non~random  Effects.  The  refractive  index  of  the  atmo- 
sphere is  a function  of  density  (hence  pressure),  temperatxire 
and  humidity.  All  of  these  parameters  have  average  values 
which  vary  in  a smooth  manner  as  functions  of  altitude  and 
position.  This  results  in  a bending  of  the  light  passing 
thzmugh  any  significantly  different  layers  of  the  atmosphere 
at  other  than  normal  incidence.  For  example,  navigators  and 
astronomers  have  long  used  a correction  for  this  refractive 
effect  for  the  viewing  of  celestial  objects  away  from  the 
zenith. 

In  principle,  much  the  same  effect  is  present  when  a 
beam  propagates  over  a path  which  includes  different 
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temperature  regions*  even  at  an  essentially  constant  alti- 
tude, For  example,  a beam  at  low  altitude  will  be  bent  when 
crossing  over  a boundary  between  woods  and  a plowed  field 
on  a sunny  day.  As  an  estimate  of  the  magnitude  of  this 
effect,  consider  a temperature  difference  of  10°  C which  is 
not  unreasonable  in  bright  sun.  An  approximate  formula 
(Ref  34il02)  for  the  index  of  refraction  in  the  optical 
region  is 


n . 1 ^ 80  X 10-^  t 
T° 


(3) 


with  T in  °K  and  p in  millibars,  A simple  application  of 
Snell's  Law,  assuming  a sharp  boundary  between  the  areas, 
yields,  for  temperatures  of  295  and  305  °K  and  pressure  of 
1000  mbar,  a change  in  angle  of  8.9  microradians  for  an 
incident  angle  of  45°,  If  this  happens  near  the  beginning 
of  a 10  kilometer  path,  the  beam  position  at  the  end  of  the 
path  will  be  moved  8,9  cm  from  its  unperturbed  position. 

For  a highly  collimated  or  focused  beam,  this  could  be  a 
significant  factor. 

For  systems  which  use  two  way  propagation,  as  with 
reflected  signals,  or  with  beacon  signals,  as  sometimes  used 
In  communications  links,  these  effects  are  not  of  importance, 
since  the  reciprocity  inherent  ih  the  geometrical  optics 
nature  of  the  effect  leads  to  equal  and  compensatory 
refraction  in  both  directions.  For  this  reason,  and  because 
the  beams  being  used  are  somewhat  divergent,  these  effects 
are  not  considered  in  this  model. 
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Random  Effects,  The  optical  refractive  index  of  the 
atmosphere  has  been  seen  to  be  a function  of  temperature  and 
pressure.  These  two  physical  parameters,  in  addition  to 
their  relatively  well-behaved  variations  described  above, 
also  have  random  components.  These  random  inhomogeneities 

t 

in  the  atmosphere  have  been  treated  in  great  detail,  in  par- 
ticular by  Tatarski  (Ref  34).  This  book  describes  the  struc- 
ture of  the  turbulent  fluctuations,  from  both  theoretical 
considerations  and  results  of  measurements.  His  treatment 
is  very  thorough  and  has  served  as  the  basis  for  a great 
deal  of  subsequent  investigation  aimed  at  verifying  or 
extending  the  treatment. 

The  basis  of  the  theory  of  the  structure  of  turbulence 
is  the  idea  that  energy  is  introduced  into  the  atmosphere  on 
a large  scale,  with  a characteristic  size  L^,  usually  called 
the  outer  scale  of  turbulence.  This  is  most  likely  some 
form  of  local  heating,  as  in  differential  heating  by  the  sun, 
depending  upon  local  topography.  This  heated  body  of  air, 
called  a turbulent  eddy,  then  breaks  up  into  smaller  eddies. 
This  continues  until  the  critical  Reynolds  number  is  reached. 
These  eddies  then  are  characterized  by  a size  the  inner 
scale  of  turbulence.  They  no  longer  break  up,  but  lose  their 
excess  energy  through  viscous  dissipation  (Ref  34i50-54). 
These  turbulent  eddies  are  the  source  of  the  random  effects 
on  propagation. 

The  results  of  applying  a small  random  perturbation  to 
the  refractive  index  (or  the  dielectric  constant)  in  the 
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solution  of  the  usual  wave  equation  lead  to  variations  in 
irradiance,  phase  and  angle  of  arrival,  A great  deal  of 
Tatarski’s  book  (Ref  34)  and  most  of  the  other  theoretical 
literature  on  turbulence  is  concerned  with  improving  this 
approach  and  solving  the  resulting  equations,  or  with 
assuming  and  testing  various  distributions  of  the  turbulent 
eddies.  The  various  effects  are  described  below  and  the 
approximations  and  formulas  used  in  the  model  are  discussed 
in  Chapter  V. 

Irradiance  Scintillations.  The  predominant  effect  of 
atmospheric  turbulence  on  a light  beam  is  that  of  irradiance 
scintillation.  This  effect  is  what  makes  stars  twinkle,  and 
has  been  a problem  to  astronomers  for  many  years.  Tatarski 
has  demonstrated  that  the  effect  can  be  considered  as  a 
multiple  scattering  phenomenon  (Ref  34*166-177).  A great 
deal  of  his  book,  and  a great  deal  of  work  by  others,  has 
gone  into  understanding  the  nature  of  the  effect, 

A method  based  upon  geometrical  optics  was  first  used 

to  try  to  understand  propagation  through  a random  medium. 

This  was  only  valid  for  very  short  paths,  on  the  order  of 
2 

ki©  , where  k is  the  wavenumber  of  the  signal.  Hence,  since 
i©  in  the  lower  atmosphezo  is  on  the  order  of  millimeters, 
this  approach  was  valid  only  for  paths  of  a few  hundred 
meters,  at  most. 

The  approach  used  by  Tatarski  was  based  upon  an  approxi 
mation  attributed  to  Rytov.  In  this  approximation,  a 

IS 
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perturbation  method  is  applied  to  the  natural  logarithm  of 
the  electric  field.  That  is,  the  electric  field  E is 
expressed  as 

E *»  exp(4f)  (^) 

and  then  is  expressed  as 

where  is  the  part  of  i)j  which  is  the  vacuum  result.  The 
assumption  that  is  a small  perturbation  is  then  made,  and 
solutions  for  can  be  obtained.  Tatarski  calls  this  the 
Method  of  Smooth  Perturbations  (Ref  34 i 218-258).  The  range 
of  validity  of  this  Rytov  approximation  is  limited  by  the 
assumption  that  « lil/gj.  It  was  found  in  the  late 

1960*8,  that  for  paths  longer  than  about  one  kilometer, 
measured  data  did  not  agree  well  with  predictions. 

Since  that  time,  many  efforts  have  been  made  to  extend 
the  theory  to  the  case  of  longer  paths  and  for  stronger  tur- 
bulence. Fante  (Ref  lOi 11-15,  34-38)  gives  a relatively 
simple  and  straightfoirward  discussion  of  some  of  the  methods 
and  approximations. 

One  of  the  most  tenacious  problems  in  the  whole  area 
of  turbulence  has  been  the  saturation  of  log-irradiance. 

The  fluctuations  in  the  irradiance  are  usually  characterized 
by  the  variance  in  the  log  of  the  Irradiance.  This  is  mostly 
a carryover  from  all  the  work  done  with  the  Rytov  approxi- 
mation. It  is  found  (Ref  17il524)  that  the  log-intensity 


14 


variance  increases  with  path  length  and  turbulence  strength 
until  it  reaches  a saturation  value.  Above  that  path  length 
or  turbulence  strength,  the  variance  no  longer  increases. 

It  actually  tends  to  decrease  slightly. 

Many  attempts  have  been  made  to  explain  the  saturation 
effect.  None  has  been  widely  accepted  by  other  researchers. 
Yura  has  recently  (Ref  40)  presented  a model  based  upon 
Tatarski's  model  with  the  addition  of  diffraction  and  the 
loss  of  lateral  spatial  coherence  as  the  beam  propagates. 

This  model  predicts,  at  least  qualitatively,  all  of  the 
observed  effects  in  the  saturation  region. 

Another  open  question  in  the  understanding  of  the  irra- 
diance  scintillations  is  the  actual  probability  distribution 
of  irradiance.  For  weak  turbulence,  it  has  been  generally 
accepted  (Ref  10j42)  that  the  probability  distribution  is 
log-normal.  This  means  that  the  log  of  the  irradiance  has  a 
normal  distribution.  It  has  also  been  noted  that  strong 
turbulence  leads  to  the  same  distribution.  In  the  intermedi- 
ate area,  however,  there  are  significant  deviations  from 
log-normality. 

It  has  been  shown  (Ref  37)  that  the  distribution  cannot 
be  everywhere  log-normal.  Other  possible  distributions  being 
proposed  include  a normal  distribution  (Ref  3),  small  pertur- 
bations on  the  log-normal  (Ref  38)  and  a three -parameter 
truncated  log-normal  (Ref  9).  The  last  of  these  has  been 
rather  arbitrarily  chosen  for  the  present  model.  Both 
distributions  are  described  in  Chapter  V. 
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Angle  of  Arrival  Variations 


When  a wave  propagates  through  a vacuum,  the  wavefront 
remains  uniform.  As  a wave  passes  through  the  atmosphere, 
however,  different  parts  of  the  wavefront  experience  dif- 
ferent phase  shifts,  resulting  in  a non-uniform  wavefront. 
The  effect  of  this  distortion  upon  the  angle -of -arrival  at 
the  receiver  (or  at  the  target)  depends  in  part  upon  the 
definition  of  angle-of -arrival  (Ref  17 i 1542).  Using  the 
physical  description  of  angle-of -arrival  as  the  direction  of 
the  normal  to  the  surface  of  constant  phase,  Fante  (Ref  10 t 
44-48)  has  defined  a as  the  local  angle  of  arrival,  the 
difference  between  the  normal  to  the  perturbed  constant 
I^se  surface  and  the  normal  to  the  unperturbed  (vacuum) 
phase  surface  (Fig.  1).  He  then  derives  an  expression  for 
the  mean- square  angle -of -arrival  flucuation. 


r*  C^(ic')dx' 

<a2>.A  


I 


(6) 


where  x is  the  source  to  receiver  range,  D*  is  the  receiver 
diameter  and  A is  a coefficient  of  order  unity,  depending 
upon  the  strength  and  scale  sizes  of  the  turbulence.  is 
the  refractive  index  structure  parameter,  discussed  on 
page  39.  This  expression  should  be  valid  for  any  non-focused 
beam. 

Using  an  effective  value  of  cj  (constant  over  the  path), 
a 10  kilometer  range,  a 10  cm  receiver  diameter,  and  A > 3 
(to  give  largest  likely  fluctuations),  we  find  the  values 
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of  V<a^>  In  Table  I. 

Table  I 

Range  of  Arrival  Angle  Variations 


(radians) 

5,0  X 10’^2 

5,7  X lO”'^ 

5,0  X 10'^^ 

• 

o 

X 

00 

• 

5,0  X 10’^^ 

5,7  X 10'^ 

5,0  X 10"^^ 

1,8  X 10"^ 

It  should  be  noted,  as  discussed  in  Chapter  V,  that 

2 

these  values  of  correspond 

to  turbulence  levels  ranging 

from  strong  to  extreme.  Even  in  the  most  extreme  case,  the 
average  (RMS)  fluctuation  is  seen  to  be  of  the  order  of 
tenths  of  milliradians.  Since  the  equipment  being  consid- 
ered has  a field-of-view  at  least  an  order  of  magnitude 
larger  than  this,  the  angle-of -arrival  fluctuation  will  be 
ignored. 

For  a more  highly  collimated  beam  and  correspondingly 
narrower  field-of-view  of  the  receiver,  angle  of  arrival 
fluctuations  could  cause  significant  variations  in  received 
power,  particularly  for  certain  combinations  of  wavelength 
and  receiver  size. 

Phase  Fluctuations,  Except  in  the  case  of  interferometric 
measurements  of  some  sort,  the  significance  of  phase 
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fluctuations  is  essentially  limited  to  the  angle-of -arrival 
effect  described  above.  Phase  fluctuations  will  therefore 
not  be  considered. 


Wander.  There  is  yet  another  random  effect  which  is 
particularly  pertinent  to  a narrow  beam,  unlike  the  other 
effects  which  affect  any  sort  of  wave.  Consider  a narrow 
beam  passing  through  a series  of  cells  of  different  sizes 
and  refractive  indices.  These  will  act  like  a series  of 
lenses.  If  the  cell  (turbulent  eddy)  is  large  compared  to  the 
beam  diameter,  the  beam  will  be  deflected,  as  in  Fig.  2.  If, 
on  the  other  hand,  the  eddy  is  small,  it  will  tend  to  broaden 
the  beam,  as  in  Fig.  3. 

In  strong  turbulence,  perhaps  because  there  are  just 
more  turbulent  eddies,  or  perhaps  because  they  are  more 
sharply  defined,  the  beam  may  even  break  up  into  multiple 
beams.  Various  studies  have  been  done  to  investigate  these 
effects  (Ref  lOt 18-24).  Even  though  a rigorous  description 
of  the  mechanism  is  not  available,  it  is  possible  to  predict 
the  sort  of  wander,  spreading  or  splitting  by  a comparison  of 
the  various  parameters.  This  is  described  in  Chapter  V, 
also. 
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Ill,  Modeling  of  Molecular  Effects 

The  molecular  attenuation  effects  in  the  atmosphere  are 
those  due  to  the  gaseous  species  present  in  the  atmosphere. 

The  gas  molecules  both  absorb  and  scatter  electromagnetic 
radiation.  The  scattering  is  essentially  classical  Rayleigh 
scattering  from  small  particles,  while  the  absorption  is  a 
quantum  mechanical  effect  which  is  strongly  wavelength 
dependent.  The  wavelengths  at  which  absorption  takes  place 
depend  upon  the  energy  levels  of  each  molecular  species. 

Hence,  the  total  absorption,  unlike  the  scattering,  is  a 
strong  function  of  the  relative  concentrations  of  the  various 
molecular  species. 

Distribution  of  Molecular  Species 

The  molecules  assumed  to  be  the  normal  constituents  of 
the  atmosphere  fall  into  two  classes.  One  of  these  is  called 
the  uniformly  mixed  gases.  These  are  nitrogen,  oxygen,  carbon 
dioxide,  carbon  monoxide,  nitrous  oxide,  and  methane.  As  the 
name  implies,  these  gases  are  assumed  to  be  uniformly  mixed 
throughout  the  atmosphere,  with  concentrations  as  shown  in 
Table  II.  These  are  the  values  used  in  the  published  reports 
of  the  Air  Force  Geophysical  Laboratory  (Refs  19,  20,  22). 
These  values  are  based  upon  extensive  measurements  over  many 
years  (Ref  22),  and  represent  average  values,  A discussion 
of  the  variability  of  these  ratios  is  found  in  Ref  23. 
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Table  II 


Mixing  Ratios  of  the  Uniformly  Mixed  Gases 


Gas 

Mixing  Ratio 

(parts  per  million,  by  volume) 

7.905  X 10^ 

^2 

2.095  X 10^ 

“2 

330 

CO 

.075 

N20 

0.28 

CH^ 

1.6 

(Ref  29 1 10) 


The  other  class  of  gases  consists  of  water  vapor  and 
ozone,  which  vary  greatly  in  concentration  with  altitude, 
goegrai^ical  location,  and  temperature.  Hence,  these  two 
gases  are  treated  separately  \dien  describing  the  atmosphere. 

Note  that  no  account  is  taken  of  other  gases  present  in 
the  atmosphere.  There  are,  first  of  all,  several  other  gases 
normally  present  in  "tiniformly  mixed"  proportions,  such  as 
the  noble  gases.  None  of  these  is  present  in  large  enough 
quantities  to  significantly  affect  absorption.  There  is  also 
the  problem  of  pollutants,  such  as  hydrocarbons.  These 
could,  particularly  in  urban  or  industrial  areas,  result  in 
significeuit  absorption  at  particular  wavelengths. 

The  very  large  number  of  possible  species  present  in 
this  sort  of  environment,  and  the  virtual  impossibility  of 
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predicting  just  what  mixture  will  be  present  in  a given  case, 
make  accounting  for  this  a virtual  impossibility.  In  any 
case,  the  sort  of  situation  in  which  these  gases  would  be 
present  in  significant  quantities  would  almost  certainly  be 
one  in  which  particulate  pollutants  from  the  same  sources 
would  provide  the  predominant  loss  mechanism.  Therefore, 
such  gases  are  generally  not  considered,  and  are  not 
included  here. 

Standard  Atmospheres 

In  order  to  provide  a common  basis  for  atmospheric 
description,  the  U.S.  Standard  Atmospheric  Supplements,  1956, 
and  the  Handbook  of  Geophysics  and  Space  Environment 
(Ref  35)  are  generally  used.  These  have  been  modified  by 
AFGL  to  include  water  vapor  and  ozone  densities  (Ref  22) . 
These  models  include  pressure,  temperature,  density,  and  the 
densities  of  water  vapor  and  ozone  as  a function  of  altitude. 
From  the  density,  the  number  densities  of  the  uniformly  mixed 
gases  can  be  calculated.  Likewise,  the  number  densities  of 
water  vapor  and  ozone  can  be  calculated  from  their  densities. 
These  number  densities  are  then  used  in  calculating  the 
effects  of  the  gases  on  radiation. 

The  six  models  commonly  used  (especially  by  AFGL)  are 
the  1962  U.S.  Standard  Atmosphere,  Tropical,  Midlatitude 
Summer  (45°  N,  July),  Midlatitude  Vi/inter  (45°  N,  January), 
Subartic  Summer  (60°  N,  July),  and  Subartic  Winter  (60°  N, 
January),  Table  III  shows  the  U.S.  Standard  Atmosphere, 
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1 


U.  S.  STANDARD  ATMOSPHERE.  1962 

Ht. 

(km) 

Pressure 

(mb) 

Temp. 

(Ok) 

DensiW 

(g/m*) 

Water  Vapor 
(g/m^) 

Ozone 

(g/vnr) 

0 

1.  013E+03 

288. 1 

1.  225E-I-03 

5.  9E+00 

5.  4E-05 

1 

8.  986E-f02 

281.6 

1.  lllE+03 

4.  2E-tOO 

5.  4E-05 

2 

7.  950E4^02 

275. 1 

1.  007E+03 

2.  9E-<-00 

5.4E-05  ' 

3 

7.  012E+02 

268.7 

9.  093E-f-02 

1.  8E-I-00 

5.  OE-05 

4 

6. 166E-I-02 

262.2 

8. 193E-I-02 

1.  lE+00 

4.  6E-05 

5 

5.  405E-(-02 

255.7 

7.  364E+02 

6.  4E-01 

4.  5E-05 

6 

4.  722E+02 

249.2 

6.  601E+02 

3.  8E-01 

4.  5E-05 

7 

4.  lllE-l-02 

242.7 

5.  900E+02 

2.  lE-01 

4.  8E-05 

8 

3.  565E-I-02 

236.2 

5.  258E-I-02 

1.  2E'01 

5.  2E-05 

9 

3.  080E-I-02 

229.7 

4.  671E+02 

4.  6E-02 

7.  lE-05 

10 

2.  650E+02 

223.2 

4.  135E-i-02 

1.  8E-02 

9.  OE-05 

11 

2.  270E-t^02 

216.8 

3.  648E-t-02 

8.  2E-03 

1.  3E-04 

12 

1.  940Et02 

216.6 

3. 119E+02 

3.  7E-03 

1.  6E-04 

13 

1,  658E+02 

216.6 

2.  666E->-02 

1.8E-03 

1.  7E-04 

14 

1.  417E+02 

216.6 

2.  279E+02 

8.  4E-04 

1.  9E-04 

15 

1.  211E-^02 

216.  6 

1.  948E+02 

7.  2E-04 

2.  lE-04 

16. 

1.  035E+02 

216.6 

1.  665E-I-02 

6.  lE-04 

2.  3E-04 

17 

8.  850E-I-01 

216.6 

1.  423E+02 

5.  2E-04 

2.  8E-04 

18 

7.  565E+01 

216.  6 

1.  216E-I-02 

4.  4E-04 

3.  2E-04 

19 

6.  467E+01 

216.  6 

1.  040E-t-02 

4.  4E-04 

3.  5E-04 

. 20 

5.  529E-t^01 

216.  6 

8.  891E-I-01 

4.  4E-04 

3.  8E-04 

21 

4.  729E^01 

217.6 

7.  572E+01 

4.  8E-04 

3.  8E-04 

22 

4.  047E-I-01 

218.6 

6.  451E-t-01 

5.  2E-04 

3.  9E-04 

23 

3.  467E+01 

219.6 

5.  500E-I-01 

5.  7E-04 

3.  8E-04 

24 

2.  972E+01 

220.6 

4.  694E-I-01 

6.  IE- 04 

3.  6E-04 

25 

2.  549E-f01 

221.6 

4.  008E-I-01 

6.  6E-04 

3.  4E-04 

30 

1. 197E+01 

226.5 

1.  841E+01 

3.  8E-04 ' 

2.  OE-04 

35 

5.  746E+00 

236.^ 

8.  463E+00 

1.  6E-04 

1.  lE-04 

40 

2.  871E+00 

250.4 

3.  996E-I-00 

6.  7E-05 

4.  9E-05 

45 

1.  491Et00 

264.2 

1.  966E400 

3.  2E-0S 

1.  7E-05 

SO 

7.  978E-01 

270.6 

1.  027E+00 

1.  2E-05 

4.  OE-06 

70 

5.  520E-02 

219.7 

8.  754E-02 

1.  5E-07 

8.  6E-08 

100 

3.  008E-04 

210.0 

4.  989E-04 

. 1.  OE-09 

4.  3E-11 

(lef  22) 

Table  III.  U.S.Steodard  Atmosphere  (1962) 


Table  IV 

Model  Atmosphere  Sea  Level  Values 


Model 

Pressure 

(mbar) 

Temp 

(0|C) 

Water -Vapor 

Ozone 

(g/m^) 

U.S.  Stand 

1013. 

288.1 

5.9 

5.4  E-05 

Tropical 

1013. 

300.0 

19.0 

5.6  E-05 

Mid-Lat  Sum 

1013. 

294.0 

14.0 

6.0  E-05 

Mid -La t Winter 

1018. 

272.2 

3.5 

6.0  E-05 

Sub-Arctic  Sum 

1010. 

287.0 

9.1 

4.9  E-05 

Sub-Arctic  Winter 

1013. 

257.1 

1.2 

4.1  E-05 

1962 ( as  an  example  of  the  models.  A comparison  of  the  sea 
level  values  is  given  in  Table  IV  as  an  indication  of  the 
readily  available  range  of  predictions.  The  water  vapor 
density  in  the  U.S.  Standard  atmosidiere  corresponds  to 
approximately  50%  relative  humidity*  while  the  other  five 
all  represent. about  80%  relative  humidity  at  sea  level. 

The  altitude  dependence  of  humidity  also  differs  (Ref  29 t 9). 
Proper  selection  of  an  atmospheric  model  is  therefore  partic- 
ularly important  %dienever  concerned  about  an  area  of  the 
spectrum  where  absorption  by  water  vapor  may  be  significant. 

Molecular  Scattering 

The  scattering*  as  stated  in  Chapter  II,  is  essentially 
Rayleigh  scatter,  and  can  be  described  by  the  formula 

a.  - (9.87  X 10'2°)  X (7) 
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where  v * ^ is  the  wavenumber  of  the  radiation.  The  con- 
stant in  this  formula  Is  an  empirical  one,  obtained  from 
LOWTRAN  3 (Ref  29ill),  This  calculation  is  performed  in  a 
subroutine  called  RAYLEE,  even  though  it  is  essentially  only 
one  step  long.  This  is  done  to  be  consistent  with  the  other 
calculations. 

Molecular  Absorption 

The  problem  of  molecular  absorption  is  much  more  diffi- 
cult to  handle.  The  Air  Force  Geophysical  Laboratory  has 
developed  several  computer  codes  designed  to  calculate  the 
molecular  absorption.  The  results  of  one  of  these  programs 
(Ref  7)  in  the  area  of  the  1.06^  wavelength,  is  shown  in 
Fig,  4,  Transmittance  is  plotted  as  a function  of  frequency 
(in  cm”^).  The  laser  line  center  frequency  is  almost  exactly 
9400.2  cm~^,  and  the  width  of  the  laser  line  is  approximately 
3 cm~^.  This  calculation  was  made  for  a midlatitude  winter 
standard  atmospheric  model  and  for  a range  of  100  km,  in 
order  to  show  the  weak  absorption  lines  present  in  the 
atmosphere  in  this  wavelength  region. 

For  both  the  1.06(i  and  .53|xm  beams,  the  absorption  coef- 
ficient is  set  equal  to  .001,  corresponding  to  a transmittance 
in  excess  of  99%  for  a 10  km  path.  This  is  based  upon  the 
calculations  from  the  AFGL  programs.  All  of  these  programs, 
as  described  in  Appendix  C,  require  significant  amounts  of 
computer  core  space  and  processor  time,  in  addition  to  lengthy 
files  of  molecular  absorption  line  parameters.  Therefore, 
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Results  of  Program  HIRACC 


this  figure  has  been  set  in  as  a conservative  measure  of 
the  absorption. 

The  absorption  coefficient  for  the  dye  laser  beam  is 
also  set  to  .001,  since  there  is  very  little  absorption  over 
most  possible  wavelengths,  and  the  calculations  necessary  to 
obtain  exact  values  are  excessive.  As  a cautionary  note, 
however,  it  should  be  understood  that  there  can  be  signifi- 
cant absorption  in  the  wavelengths  between  the  visible  and 
the  1.06  micron  beam.  In  particular,  there  is  significant 
absorption  around  .70  microns,  .83  microns,  and  between 
.89  and  .99  microns  (Ref  20).  As  new  dyes  become  available, 
these  absorption  problems  must  be  evaluated.  The  program  is 
set  up  to  flag  any  wavelength  in  these  regions,  warning  that 
more  detailed  calculations  of  absorption  are  required. 

A factor  vdiich  may  be  of  importance  for  laser  beams  of 
very  narrow  optical  bandwidth,  when  calculating  molecular 
absorption,  and  which  can  generally  be  ignored  for  any  of  the 
other  calculations,  is  that  of  precisely  defining  the  wave- 
length. All  the  absorption  compilations  are  based  upon 
wavenumbers,  which  is  essentially  a frequency  and  hence 
independent  of  medium.  The  reciprocal  of  wavenumber,  however, 
gives  the  vacuum  wavelength  of  the  radiation.  In  contrast, 
laser  manufacturers  generally  specify  the  wavelength  as 
determined  by  laboratory  measurements  of  wavelength  in  air. 
Hence,  the  reciprocal  of  1,0638  microns  would  lead  to  a 
wavenumber  of  9400.26  cm~^.  If  the  actual  vacuum  wavelength 
is  (1.0638)  X (1.0004)  * 1,0642  microns,  the  corresponding 


actual  wavenumber  Is  9396.50  cm"^.  This  Is  more  than  enough 
difference  to  move  from  the  far  wings  of  an  absorption  line 
to  the  peak.  In  the  case  of  a laser  such  as  the  iVdiYAG,  the 
llnewidth  of  the  transition  is  on  the  order  of  3 cm“^  at 
room  temperature  (Ref  32*362),  so  that  the  presence  of  one 
or  two  small  narrow  absorption  lines  may  safely  be  ignored 
and  there  is  no  need  to  be  concerned  with  the  wavelength 
ambiguity. 


IV.  Model ins  of  Aerosol  Effects 


As  has  been  stated  above,  it  is,  in  principle,  possible 
to  accurately  calculate  the  total  aerosol  attenuation  for 
the  case  of  single  scatter,  given  the  distribution  of  aero- 
sol particle  sizes  and  complex  refractive  indices.  Some  of 
the  recent  work  on  determining  and  modeling  aerosol  distri- 
butions will  be  discussed  and  the  simplifying  approximations 
used  in  the  present  model  will  be  described. 

Atmospheric  Models 

Despite  considerable  effort  by  many  researchers  over 
the  past  several  decades,  there  is  still  no  reliable  method 
of  predicting  the  particular  distribution  of  aerosol  concen- 
tration to  be  found  on  a given  day  (Ref  31).  Indeed,  there 
is  no  way  of  knowing  this  distribution  even  when  given  a full 
set  of  meteorological  data  as  normally  reported.  Many  dif- 
ferent aerosol  distributions  could  lead  to  the  same  reported 
visibilities. 

This  lack  of  precise  information  concerning  the  particle 
distribution  has  led  to  many  attempts  to  develop  models  of 
the  "typical"  atmosphere  (Refs  28,  29,  31). 

The  Air  Force  Cambridge  Research  Laboratories  (AFCRL) 
(now  renamed  the  Air  Force  Geophysical  Laboratory)  has  con- 
tinuously worked  on  improving  and  updating  these  aerosol 
models.  The  latest  generation  of  these  models  was  published 
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in  1976  (Ref  31).  This  series  of  models  consists  of  three 
boundary  layer  (up  to  2 kilometer  altitude)  models,  two 
upper  troposphere  models,  four  stratosphere  models,  and  two 
upper  atmosphere  models. 

For  the  case  of  near-earth,  nearly  horizontal  paths, 
only  the  boundary  layer  models  are  considered.  These  three 
models  are  the  rural,  urban,  and  maritime  models.  The  rural 
model  represents  the  aerosol  conditions  found  in  continental 
areas  not  influenced  by  urban  or  industrial  areas.  The 
aerosol  is  assumed  to  be  composed  of  70%  water  soluble  and 
30%  dustlike  particles. 

The  urban  model  differs  in  that  sootlike  particles  are 
assumed  to  be  added  to  the  rural  mixture,  in  the  form  mostly 
of  combustion  products.  These  carbonaceous  particles  are 
assumed  to  comprise  35%  of  the  distribution. 

The  maritime  model  includes  salt  particles  which  result 
from  sea  spray  evaporation.  This  model  is  incomplete  in 
that  It  does  not  take  windspeed  and  relative  humidity  into 
account.  These  variables  can  be  important  in  the  first 
100  meters  or  so  above  the  surface. 

Using  these  three  models  for  the  particle  distributions, 
Mie*s  scattering  theory  was  used  to  determine  the  details  of 
the  scattering  and  absorption.  For  the  present  purposes, 
only  the  behavior  of  the  overall  aerosol  extinction  coeffi- 
cient is  of  interest.  Fig,  5 shows  the  aerosol  extinction 
coefficient  as  a function  of  wavelength. 


Rural 

0.077 

1.21 

Urban 

0,088 

0.98 

Maritime 

0.141 

0.196 

These  approximations  are  used  in  the  subroutine  AROSOL 
for  the  aerosol  extinction  coefficients  corresponding  to  a 
visibility  of  23  kilometers.  These  expressions  give 
extinction  coefficients  within  2%  of  those  calculated  from 
the  exact  models,  as  incorporated  in  LOWTRAN  3B  (Ref  30). 

The  interpolation  scheme  used  for  other  visibilities  is 
that  found  in  LOWTRAN  3 (Ref  29il4).  There  it  is  found  that 

N(2)  - - b(2)  ( 9 ) 
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WAVELENGTH  (microns) 


(Ref  31) 

Fig.  6.  Aerosol  Extinction  Approximation  (Rural  Model) 
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(Ref  31) 

Fig.  7.  Aerosol  Extinction  Approximation  (Urban  Model) 
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(Ref  31) 

Fig.  8.  Aerosol  Extinction  Approximation 
(Maritime  Model) 
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where  N(z)  is  the  aerosol  number  density  which  is  a function 
of  altitude  and  is  proportional  to  the  aerosol  extinction 
coefficient  o_. 

The  procedure  used  for  this  interpolation  in  LOWTRAN  is 
based  upon  a tabulation  of  N at  1 kilometer  altitude  inter- 
vals, for  visibilities  of  5 and  23  kilometers.  For  the 
altitudes  of  interest  a and  b are  calculated  from  £q  (9) 
above  and  then  the  value  of  o_  is  calculated.  For  the 

cl 

present  model,  only  the  ground  level  number  densities  are 
used,  yielding 

a,  - X (^■^894  X 10^  . 2 085  ^ jo!;  ^o) 

*VIS  2.83  X lO-’  \ VIS  ' 

It  should  be  noted  that  both  in  LOWTRAN  and  in  this 
analysis  the  number  densities  used  are  those  for  the  "average 
continental  aerosol  model,"  which  was  a preliminary  version 
of  the  rural  model  used  in  LOWTRAN  3B  (Ref  30).  This  should 
introduce  no  significant  errors  into  the  calculation,  how- 
ever. 

The  AROSOL  subroutine  also  includes  provision  for  a non- 

) 

standard  model  of  the  aerosol  distribution.  This  is  accom- 
plished by  entering  new  values  for  C and  D,  the  parameters 
as  fotind  in  Table  V.  This,  of  course,  assumes  that  the  new 
distribution  can  also  be  approximated  by  a straight  line  on 
the  log-log  plot  as  in  Fig.  5.  Note  that  this  would  allow 
the  subroutine  to  be  used  in  other  wavelength  intervals, 
anomalous  refractive  indices  around  6-10  microns  wavelength 
are  avoided. 
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V.  Modeling  of  Turbulent  Effects 

Several  very  basic  assumptions  are  made  concerning  the 
nature  of  the  turbulent  effects.  The  first  is  that  the 
interaction  of  the  beam  with  the  turbulent  eddies  falls  into 
the  inertial  subrange  of  turbulence.  This  requires  that  the 
Reynolds  nxjmber  be  high  enough  so  that  there  is  a significant 
gap  between  the  size  scales  where  energy  is  introduced  into 
the  atmosphere  and  where  it  is  viscously  dissipated.  Tatar- 
skii (Ref  34i51-59)  discusses  this  requirement,  along  with 
the  detailed  description  of  the  scale  lengths.  Lawrence  and 
Strohbehn  (Ref  17*1526-1527)  discuss  the  turbulence  parameters 
as  well.  From  their  discussion,  the  inner  and  outer  scale 
lengths  have  been  chosen  as  1 mm  and  50  meters,  respectively. 
Hence,  all  the  significant  turbulent  eddies  are  assumed  to 
have  sizes  between  1 mm  and  50  meters. 

The  laser  beam  is  assumed  to  have  the  initial  form  of  a 
Gaussian  beam. 

Eq(0,  p)  ■ a expC-  (^  + ^)(p^)]  (11) 

where  D is  the  diameter  of  the  beam  and  F is  the  radius  of 
curvature  of  the  wavefront. 

After  reflection  from  the  target,  the  beam  is  assumed 
to  have  the  same  form.  This  has  been  chosen  as  generally 
reasonable  for  any  size  target,  even  though  reflection  from 
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k small  comer  cube  may  more  nearly  resemble  a uniform  small 
source*  which  could  be  treated  as  a point  source  of  spherical 
waves* 

At  the  target  plane  (and  receiver  plane  for  the  reflected 
beam)  the  beam  is  assumed  to  have  an  instantaneous  Gaussieui 
shape*  with  size  parameters  discussed  below  under  beam  wander. 


Turbulence  Strength 

2 

The  parameter  C^*  called  the  refractive -index  structure 

2 

parameter*  is  the  measure  of  the  strength  of  turbulence. 

-k  • — ' 

has  units  of  meters  It  is  actually  a measure  of  the 

difference  in  refractive-index  between  two  nearby  points. 
Actual  measurements  of  the  parameter  are  always  done  indi- 
rectly* at  least  in  the  optical  wavelength  range.  It  is 

either  done  by  observing  the  optical  turbulence  effects  and 

2 2 
going  back  to  calculate  C^*  or  by  taking  measurements  of  C^* 

the  equivalent  parameter  for  temperature  variations.  The 

2 2 
value  of  can  then  be  simply  related  to  by  the  formula 

(Ref  34 1 102) 


80  X 10 

ZT 


-6 


(12) 


«diere  the  numerical  coefficient  of  80  is  actually  slightly 

wavelength  dependent. 

2 

In  general*  is  strongly  dependent  upon  position  and 

2 

altitude.  Various  expressions  have  been  used  to  express 
as  a function  of  altitude  (actually  height  above  ground).  A 
comon  model  is  that  of  Hufnagel  (Ref  13)  which  is  presented 
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in  graphical  form  (Fig.  9)  rather  than  as  an  analytic  model. 
It  can  be  readily  seen  that  the  weather  conditions  strongly 
affect  the  value,  as  do  time  of  day  and  altitude.  Rather 
than  using  a particular  analytic  function  or  tabulating  data 

for  the  curves  of  the  Hufnagel  model,  an  "effective  value" 

2 2 
of  is  used  in  this  model.  A value  of  C.^  is  selected  as 

representative  of  the  path  as  a whole  and  treated  as  a con> 

stant  in  the  calculations.  The  default  value  is  set  to 

5.0  X 10“^^,  which  can  be  seen  to  correspond  to  a very  low 

value  at  the  surface  but  a high  value  at  a height  of 

one  hundred  meters,  hence  not  unreasonable  for  the  entire 

path.  Since  this  choice  is  clearly  arbitrary,  it  can  be 

varied  to  fit  the  results  of  measurements. 

Beam  Wander 

The  irfienomenon  of  beam  wander  is  treated  as  described 
Fante  (Ref  10 i 18-24).  The  beam  is  assumed  to  have  an 
instantaneous  Gaussian  cross-section  at  the  target  (or 
receiver)  plane.  The  short  term  beam  radius  is  defined 
as  the  radius  at  which  the  short  term  average  irradiance  is 
reduced  by  a factor  of  e~^  fzmm  the  maximum.  In  general, 
then,  this  short  term  spot  will  be  deflected  by  the  large 
eddies  in  the  beam.  The  long  tern  average  irradiance  will 
therefore  be  spread  over  a larger  area.  The  e*^  point  of 
this  distribution  is  defined  as  p^.  Another  parameter  p^  is 
defined  as  the  measure  of  the  deflection,  or  beam  centroid 
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wander*  where 


2 2 2 
Pc  “ Pl  - Ps 


(13) 


Then*  using  the  waveform  of  Eq  (f/  )*  it  can  be  shown 
(Ref  10 I 21)  that 


/ 2\  ~ 4x^  X 

\pl/“  +2r'^  - - 


2 4x2 

+ -y-j 

^Po 


(14) 


for  X « (k2  c2  jtjiyl 


and 


k^D‘ 


1 i 

1 * rt 


2 ^ 2.2  x^ 


TJ 


(15) 


.22 

for  X » C“ 


where*  as  usual*  k is  the  wavenumber  2it/X  and  is  defined 


as 


- (.626  k2  X c2)  ^ 


(16) 


The  first  two  terms  in  each  of  the  above  expressions 
make  up  the  vacuum  beam  spread  (VBS  in  the  program)  and  the 
third  term  is  the  turbulent  contribution. 

The  calculation  of  and  is  not  so  straightforward. 
There  are  four  different  combinations  of  length  parameters 
which  lead  to  different  expressions  for  the  short  term 
effects.  These  depend  upon  the  coherence  length  p^*  defined 
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above}  beam  diameter  Di  outer  scale  of  turbulence  L^i  range 

X}  wavenumber  k and  L,  where  L Is  the  smaller  of  and  D. 

It  is  also  assumed  that  D > i . 

o 

Case  1.  Pq  ^ then 


VBS  + 


4x^ 

~T~1 


1 


(17) 


and 


^ 2,97  X 


(18) 


Case  2. 


2 

If  Pq  M D and  x < kL  , there  are  no  simple 

expressions  for  p„  and  p_.  Fante  (Ref  10 i 

s c 

23-24)  gives  a graphical  solution  in  terms 
of  a parameter  3,  where 


(19) 


The  grai^lcal  solution  leads  to  a value  of 
where 


(20) 


In  the  computer  model t the  graphical  solu- 
tion has  been  approximated  by  storing  values 
of  the  three  plotted  curves  (Fig.  10)  in  a 
two  dimensional  array  (RATIO)  and  using  an 
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(Ref  10i24) 

Fig.  10.  Parameters  for  Beam  Wander  Calculation 


Interpolation  scheme  to  find  (i.  Then 

s 

is  calculated  from  Eq  (17)*  and  from 

w 

Eq  (13). 

2 

Case  3.  If  p^  » D and  x < kL  * beam  wander  is 

negligible  and  p^  « p^.  This  corresponds 
to  short  ranges  (a  few  hundred  meters)  or 
weak  turbulence. 

2 

Case  4.  If  x » kL  * the  beam  breaks  up  into  mul- 
tiple patches.  The  instantaneous  pattern 
consists  of  several  spots*  and  the  beam 
centroid  wander  is  negligible.  An 
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approximate  expression  for  (Ref  10i23) 

2 

is  available  for  x » kD  and  D » p^. 

(pc)  "*  ^ 

This  expression  has  been  used  here  for  p^ 

2 

in  all  cases  for  x > kL  . In  this  case* 
p^  has  no  well  defined  meaning.  It  has 
been  arbitrarily  set  equal  to  VvBS,  which 
is  the  vacuum  beam  spread.  Clearly*  the 
numbers  resulting  from  calculations  using 
this  case  are  suspect.  The  calculations 
above  give  no  idea  as  to  the  number  of  beams 
or  their  locations*  only  that  the  root  mean 
square  radius  of  the  area  where  they  form  is 
Pj^.  A rather  arbitrary  cutoff  of  Pj^  > 1.5  p^ 
was  set  in  the  program  to  warn  that  beam 
breakup  is  causing  serious  spreading  of  the 
long  term  average.  This  means  that  the 
short  term  prediction  of  irradiance  on  axis 
is  questionable*  and  is  probably  more  vari- 
able on  a pulse-to-pulse  basis  them  the 
calculated  results  indicate. 

These  four  cases  are  considered  and  p^*  p^*  and  p^  are  cal- 
culated in  subroutine  AXIAL.  The  mean  irradiance  and  char- 
acteristic time  of  beam  wander  (WNDTIM)  are  also  calculated 
in  AXIAL. 
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Temporal  Behavior  of  Beam  Wander.  The  entire  treat- 
ment of  turbulence  has  been  based  upon  an  assumption  of  the 
"Taylor  frozen  flow  hypothesis,"  This  is  the  assumption  that 
the  turbulence  encountered  by  the  beam  is  in  the  form  of 
small  cells  which  do  not  change  while  they  are  seen  by  the 
beam.  Rather,  the  changes  in  turbulence  are  due  to  the  flow 
of  such  cells  across  the  beam  (Ref  34t88-91).  Clearly,  then, 
for  the  case  of  beam  wander  being  considered,  the  deflection 
or  spreading  of  the  beam  can  be  expected  to  change  on  a time 
scale  of  D/V*  where  D is  the  beam  diameter  and  V is  the 
transverse  wind  velocity.  Neither  of  these  parameters  is  a 
constant  for  the  case  of  a diverging  beam  in  the  real  atmo- 
sphere. An  average  value  of  V is  one  of  the  input  parameters. 
The  value  of  D used  is  the  initial  beam  diameter  at  the 
transmitter  for  the  forward  beam  and  target  diameter  for  the 
reflected  beam.  This  gives  the  shortest  time  scale  of 
changes  to  be  expected. 

It  is  to  be  noted  that  for  reasonable  values  of  D and 

V,  these  characteristic  times  (Table  VI)  are  of  the  same 

order  of  magnitude  or  less  than  the  pulse  repetition  time 

for  a low  duty  cycle  pulsed  laser  (3  x 10  seconds).  This 

means  that  the  wander  effect  can  be  expected  to  change  from 

one  pulse  to  the  next.  The  times  are  long  coapared  to  the 

-8 

length  of  a single  pulse  (10  sectxids),  however,  so  that 
each  sees  essentially  a single  set  of  cells.  Thus,  there  is 
no  time  averaging  within  each  pulse,  and  successive  pulses 
see  uncorrelated  turbulence  cells. 
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Table  VI 

Characteristic  Time  Scale  of  Beam  Wander 


D 

(m) 

V 

(m/sec) 

(v) 

(knots) 

Wander  Time  Scale 
(seconds) 

.01 

1 

(1.9) 

1 X 10-2 

.01 

2 

(3.9) 

5 X 10”^ 

.01 

3 

(5.8) 

3 X 10’^ 

.01 

5 

(9.7) 

2 X 10"^ 

.10 

2 

(3.9) 

5 X 10’2 

Mean  Irradiance 

The  calculation  of  the  mean  irradiance  on  axis  at  a 
distance  x from  the  transmitter  (or  reflector)  starting  with 
the  basic  electromagnetic  wave  equation  is  possible  using 
Fourier  transform  techniques  (Ref  10i24-25).  This  involves 
assuming  a form  for  the  wavenumber  spectrum  of  the  index  of 
refraction  fluctuations  and  for  the  input  waveform.  This 
wavenumber  spectrum  is  discussed  below  since  it  is  required 
when  considering  the  irradiance  scintillations.  For  the 
present  calculation,  a simpler  approach  is  taken. 

A Gaussian  input  beam  is  assumed.  In  terms  of  the 
wander  description  above,  it  is  assumed  that,  since  all 
refractive -index  differences  are  small,  so  are  the  deflections. 
In  terms  of  the  more  mathematical  treatment  (Ref  34)  which 
treats  all  of  turbulence  as  a scattering  from  random  inhomo- 
geneities, the  equivalent  assumption  is  that,  since  » X, 
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essentially  all  the  scattering  is  in  the  foivard  direction. 

In  either  case*  the  result  of  the  assumption  is  that 
all  of  the  transmitted  power  (after  molecular  and  aerosol 
losses  are  accovinted  for)  reaches  the  target  plane.  The 
spatial  distribution  is  assuned  to  be  Gaussieui.  The  value 
of  calculated  above  Is  tciken  as  the  e'^  radius  for  the 
distribution.  Hence*  if  A is  the  peak  irradiance,  the 
distribution  is  assumed  to  be 

Eg(p)  - A exp[-  (p/p^)^]  (22) 

and  the  total  power  in  the  target  plane  is 

'’l  ‘ toil  " Ee(p)  ^ (23) 

area 

Pj  ■ 2itA  r p expC-  (p/pj)^]  dp  (24) 

*^o 

- IT  A p2  (25) 

Using  the  conservation  of  energy  argument  as  stated 
above*  this  is  set  equal  to  the  input  power*  after  atten- 
uation* and  solved  for  A*  the  average  irradiance  on  axis 
(called  AXINT  in  the  program). 

Irradiance  Scintillations 

By  far  the  largest  effect  of  turbulence  on  these 
diverging  beams  is  the  scintillation  effect.  Any  beam  trhich 
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passes  through  more  than  a few  meters  of  turbulent  atmosphere 
suffers  redistribution  of  its  energy.  The  observed  irradiance 
then  displays  fluctuations  which  are  called  scintillations. 

The  phenomenon  has  been  studied  both  experimentally  and 
theoretically.  Tatarski's  woirk  (Ref  34)  is  the  standard 
reference  on  the  basic  theory.  He  also  includes  a great  deal 
of  comparison  of  theory  with  the  experimental  data  which  was 
available. 

The  theoretical  treatment  (Ref  34 i 166-290)  is  essentially 
complete  for  plane  and  spherical  waves,  although  approxi- 
mations are  necessary  in  order  to  solve  the  equations.  It  is 
first  asstimed  that  the  scalar  wave  equation  can  be  used, 
rather  than  the  full  vector  equations.  This  requires  that 
polarization  effects  be  insignificant.  This  is  a quite  good 
• asstjmption  for  X « The  wave  eqviation  then  takes  the 

form 

7 ^ + k^(l  + ep  ijf  • 0 (26) 

where  i|r  is  the  field,  k the  wavenumber  and  the  small  ran- 
dom part  of  the  permittivity. 

The  scalar  wave  equation  is  then  solved  using  the  method 
of  smooth  perturbations,  and  using  what  is  generally  called 
the  Rytov  approximation.  In  this  case,  the  scalar  wave 
equation  is  re%nritten  in  terms  of  the  log  of  i|r. 

7^  * + (7  ®)^  + k^  + k^€j  ■ 0 (27) 
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This  non-linear  equation  cannot  be  solved  exactly  and  con- 
tains the  random  function  The  equation  is  solved  by  a 

perturbation  method  in  powers  of  v.  where  v ■>  ■wf)- 

After  considerable  effort,  it  is  shown  that  the  result 
ceui  be  considered  as  a case  of  multiple  scattering  from  the 
random  inhomogeneities  of  the  medium.  After  making  several 
approximations  which  limit  the  range  of  applicability  of  the 
results,  expressions  for  the  variance  in  the  log-amplitude 
are  obtained.  The  log-amplitude  is  related  to  $ by  the 
expression  $ ■ X + 1-S  ) i.e. , Re($).  The  expression 

for  the  so  called  log-variance  arei 


Plane  waves 

- .308  (28) 

Spherical  waves 

- .124  k^  x^  (29) 

The  conditions  for  validity  of  the  above  expressions  includes 

A X » i© 

x*^  L©"^  » 1 
and 


X « i©^ 


Tatarskii  then  states  (Ref  34s 258)  that  experimental  evidence 
indicates  that  his  results,  as  above,  are  Justified  for 
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Many  workers  have  continued  the  theoretical  work  and 

experiment  verification.  Lawrence  and  Strohbehn  (Ref  17) 

summarized  the  status  of  investigation  into  the  probability 

distribution  of  irradiance  and  the  calculation  of  the  log- 

variance  in  1970,  They  tabulate  the  formulas  available  for 

calculation  of  the  log-variance  under  various  assumptions. 

Listed  below  are  only  those  pertaining  to  this  study.  A 

homogeneous  medium  is  assumed,  calling  for  the  choice  of  an 

2 

effective  value  of  C^.  The  Kolmogorov  spectrum  is  again 

2 

assumed.  Ranges  are  assumed  to  be  such  that  x » which 

is  not  very  restrictive.  The  beam  wave  formula  is  from 
Ishimaru  (Ref  14).  The  equations  here  are  for  the  log- 
variance  of  irradiance  (variance  in  the  log  of  irradiance), 
rather  than  amplitude.  These  are  simply  related  by  t.he  fact 
that  the  izxadiance  is  Just  the  square  of  the  amplitude,  so 
that 


*ln. 


4 ot 


(30) 


This  is  the  LOGVAR  used  in  the  program,  and  is 

^ 2 

hereafter  called  a • 


Plane  wavet 


Spherical  wavei 


- 1.23 


^ - 0.50  Cj  ^ 


(31) 


(32) 


51 


Gaussian  beam  waves 


,2  . 4„2^q  Q33  c|)r(-  5/6)k^ 


a^x 


1 + a^x 


^iFi  -%.ll  ^ 


Mil 


w^[l+(a^x)‘‘] 


2-  - gCa^x) 


(33) 


where  w » Initial  beam  radian  (DlAM/2) 
o 

2 

» X/(tt  Wq) 

(a,ctz)  is  the  Kummer  function 
gCttjX)  » Re 


^l+io^x 


HI  ia^x 


2^1 1 'I' 


2F|^(a,bictz)  is  the  hypergeometric  function. 

As  ctirrently  set  up»  the  program  calculates  the  log -variance 
for  the  Gaussian  beam  wave  case.  The  distance  p is  the 
distance  off  the  axis  of  the  beam  at  which  the  log -variance 
is  calculated.  Only  the  on-axis  case  is  considered,  so 
that  p * 0. 

Experimental  results  have  indicated  that  the  log-variance 
saturates  at  a value  of  about  2.4,  and  indeed  decreases 
slightly  (Ref  17 i 1524).  Fante  (Ref  10 i 34-35)  discusses 
briefly,  and  gives  extensive  references  to,  the  large  volume 
of  effort  which  has  gone  into  explaining  the  behavior  in 
strong  turbulence.  He  then  presents  an  approximate  formula 
for  the  log- variance  in  strong  turbulence. 


A 


2 0 
where  is  the  nomallzed  variance  in  the  lrradl^ulce,  Is 

the  plane  wave  solution  above  (Eq  31 )»  and  strong  turbulence 

Is  defined  as  existing  when  » 1.  It  Is  easy  to  show 

that 

- ln(l  oj)  (35) 

This  gives  a plane  wave  solution  for  strong  turbulence.  The 

effects  of  scintillation  are  stronger  on  a plane  wave  rhAn  on 

a beam  or  spherical  wave,  so  that  the  calculated  log-varlauice 

will  be  stronger  than  is  appropriate.  This  should  be  a 

relatively  minor  effect  and  will  in  some  measure  compensate 

for  the  large  variation  In  Irradlance  due  to  beam  breakup 

In  strong  turbulence  which  Is  not  adequately  accounted  for. 

The  empirical  formula  is  used  for  > 1,8,  which  Is  approx- 
2 2 

Imately  where  a£  = in  the  formula. 

Aperture  Averaging 

2 

The  fluctuations  a calculated  previously  apply  strictly 

speaking  only  to  a point  receiver.  Any  finite  receiver  will 

see,  not  o^,  but  an  average  of  the  fluctuations  over  the 

aperture.  The  appropriate  average  Involves  Integrating  the 

amplitude  (or  Irradlance)  covariance  functions,  which  have 

not  been  Included  here,  over  the  aperture.  Fante  presents  a 

solution  In  graphical  form  (Fig.  11),  In  terms  of  aperture 

diameter  (Ref  10 i 38*40) , The  function  G Is  the  ratio  of 

observed  fluctuations  to  those  for  a point  aperture. 

In  the  program,  G (called  APAV)  has  been  approximated 
by  a fxinctlon  which  approximates  the  two  curves  of  Fig.  11 
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Fig,  11.  Aperture  Averaging  Approximations 


quite  well  and  is  assumed  to  be  better  than  an  interpolation 
for  an  estimate  in  the  intermediate  region. 


f r / 

).09  + 0.91  exp<-  1.0  7^  y 1 + df 


(36) 


Temporal  Behavior  of  Scintillations 

The  scintillation  behavior*  like  the  wander  effect*  has 
a frequency  spectrum  which  is  a highly  complex  fvinction  of 
the  irradiance.  Fante  (Ref  10i40-41)  refers  to  several 
approximate  solutions  for  the  spectrum  in  various  limits. 

It  has  been  experimentally  foxjnd  that  the  %ridth  of  the  fre- 
quency spectrum  can  be  expressed  (at  least  in  order  of 
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i&agnltude)  by  simple  expressions.  The  Inverse  of  the  width 
is  taken  as  the  characteristic  time  of  scintillations 
(SCNTIM). 


2 

for  < 1 

SCNTIM  = 

(37) 

2 

for  0^  > 1 

Pd 

SCNTIM  » ^ 

(38) 

.36  (Xx)^ 

"P  - (,2)V 

is  a transverse  correlation 

distance 

for  the  case  of  strong  turbulence  (Ref  I0i37). 

Probability  Distribution  of  Irradiance 

It  is  found  theoretically,  as  a direct  result  of  using 
the  Rytov  method  of  dealing  with  log  of  the  amplitude  as  the 
ftindamental  variable,  that  the  irradiance  probability  distri- 
bution is  log-normal.  This  means  that  the  log  of  the  irra- 
diance is  normally  distributed.  Experimental  results  have 
shown  good  agreement  with  this  distribution,  within  the 
limits  of  applicability  of  the  Rytov  method.  The  limit 
described  by  Fante  (Ref  I0i42)  is  that  log-normality  holds 
for  o£  < 0.3.  Purely  on  an  empirical  basis,  it  has  been 
observed  that  the  log-normal  distribution  is  also  a good 
approximation  in  the  realm  of  > 25. 

In  the  intermediate  regime,  the  expected  probability 
distribution  is  by  no  means  clear.  Fante  (Ref  10i43) 
discusses  some  of  the  distributions  proposed.  These  include 
perturbed  log-normal,  Rice-Nakagami,  Rayleigh  and  three- 
parameter  truncated  log-normal.  It  has  recently  been 
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proposed  (Ref  3)  that  a normal  distribution  is  the  most 
nearly  correct.  A basic  problem  in  experimental  verification 
is  the  wide  spread  generally  observed  in  the  data.  Many 
different  distributions  can  usually  be  made  to  fit  the  data. 

The  distribution  chosen  here  is  the  three -parameter 
truncated  log-normal  (Appendix  D).  following  Davidson  and 
Gonzales-del -Valle  (Ref  9).  This  distribution  gives  an 
excellent  fit  to  their  data.  The  measurements  were  done  in 
water  in  a carefully  controlled  experiment  in  the  laboratory. 
They  give  no  indication  of  the  appropriate  values  of  the 
offset  parameter  for  use  in  the  atmosphere,  so  it  is  arbi- 
trarily picked  as  a small  percentage  (10~^)  of  the  mean 
irradiance.  This  should  give  reasonable  results,  since  the 
observed  distributions  have  generally  been  not  too  different 
from  the  two-parameter  log-normal.  The  truncation  point  is 
clearly  at  zero  irradiance,  since  no  negative  irradiance  can 
be  observed.  This  does  not  enter  explicitly  into  the  calcu- 
lations, since  the  mean  and  variance  are  calculated  sepa- 
rately. 

The  calculation  of  the  irradiance  probability  distri- 
bution at  the  receiver  is  not  Just  a matter  of  calculating 
the  log-variance  and  mean  of  the  irradiance.  The  values  of 
log-variance  and  mean  irradiance  on  axis  at  the  target  plane 
are  first  calculated.  The  beam  wander  parameters  are  then 
applied  to  determine  the  range  of  reflected  power. 

First,  the  irradiances  on  axis  at  the  target  corre- 
sponding to  adding  or  subtracting  the  log-variance  to  the 
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log  of  the  mean  Irradiance  are  calculated,  in  subroutine 

2 

TURB.  For  the  strong  turbulence  case,  where  Is  the  cal- 
culated quantity,  it  is  applied  directly  to  the  irradiance. 
These  values  of  irradiance  are  called  UPRLIM  and  BTMLIM 
respectively. 

In  subroutine  TARGET,  the  spatial  distribution  is 
assumed  to  be  Gaussian  in  both  the  long  and  short  term 
cases,  with  axial  distance  parameters  and  p^,  as  above. 
Factors  (IL  and  TS)  are  computed  to  give  the  fraction  of 
total  power  in  the  target  plane  which  is  reflected  by  a 
target  (perfectly  reflecting)  of  diameter  TD,  for  the  given 
distribution,  with  the  pattern  centered  on  axis.  These 
reflected  powers  are  then  calculated  for  the  incident  powers 
corresponding  to  irradiances  AXINT,  UPRLIM  and  BIMLIM, 
giving  three  values  called  LTAP  (long  term  average  power) , 
UPL  and  BTL  for  the  long  term  case,  and  MOST,  UPS  and  BTS 
for  the  short  term  case.  In  the  long  term  case,  these 
represent  the  expected  range  of  power  intercepted  by  the 
target. 

In  the  short  term  case,  the  expected  wander  of  the  beam 
has  to  be  taken  into  account.  This  is  done  by  treating  p., 
the  measure  of  beam  centroid  wander,  as  the  average  instan- 
taneous, off-axis  position  of  the  beam  center.  Applying 
this  to  the  calculated  on-axis  values,  and  assuming  a small 
target,  gives  a factor  of  exp[-  (p„/p_)^]  to  be  applied  to 
MOST  and  BTS  to  give  the  average  and  lower  expected  powers, 
STAP  and  BTSC.  The  upper  power  is  not  reduced,  in  order  to 
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include  the  values  of  power  corresponding  to  little  wander 
in  the  distribution.  The  reflectivity  is  applied  to  all 
the  final  powers  to  give  reflected  power  values,  (These  are 
named  with  an  R in  the  program;  e.g. , BTS,  when  reflectivity 

is  considered,  becomes  BTSR.)  i 

] 
5 

A logical  sorting  depending  upon  the  times  character- 

] 

istic  of  the  wander,  the  scintillation  and  the  pulsed  beam  | 

is  then  performed  to  determine  the  proper  values  of  expected 
power  to  select.  Fig.  12  shows  the  options.  AP  is  the 
average  reflected  power,  BP  is  the  lower  variance  point 
power  and  UP  is  the  upper  variance  point  power.  In  most 
cases  of  interest,  the  wander  and  scintillation  times  both 
fall  between  the  pulse  duration  and  the  pulse  repetition 
time.  The  powers  used  are  then  LTARP,  UPSR  and  whichever  is 
lower  between  BTSCR  and  BTLR. 

The  three  selected  values  of  power  are  sent  to  sub- 
routine RTNPTH.  There,  they  are  used  in  three  calls  to  sub- 
routine TURB  (after  transmittance  losses  are  applied).  Each 
of  these  calls  to  TURB  returns  an  average  value  of  axial 
irradiance.  These  are  referred  to  as  ARPP,  URPP  and  BRPP. 

The  differences  (VRPP--ARPP)  and  (ARPP- -BRPP)  are  evaluated 
and  the  larger  of  the  two  is  taken  as  a variance  measure, 
from  which  a log- variance  measure  VW  is  calculated  using 
Eq  (35). 

This  log-variance  VW  is  associated  with  both  beam  wander 
and  scintillations.  When  wander  is  not  important,  it  need 
not  be  considered.  Therefore,  throughout  the  wander 
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Fig.  12.  Relationship  of  Characteristic  Times 


Fig.  12.  (continued) 
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calculatlonst  whenever  the  wander  is  found  to  be  insignifi- 
cant • either  because  it  is  time  averaged  or  because  the 
magnitude  is  small,  an  index  JERKY  is  set  equal  to  zero. 

When  JERKY  is  zero  in  the  probability  distribution  calcu- 
lations, VW  is  not  used.  If  JERKY  is  not  equal  to  zero,  the 
larger  of  VW  and  the  log-variance  firom  the  scintillation 
calculation  i s used. 

The  log-variance  due  to  scintillations  alone  is  calcu- 
lated for  the  receiver  location  simply  by  adding  the  vari- 
ances from  both  paths.  The  variance  used  for  the  return 
path  is  the  one  calculated  in  TURB  along  with  the  average 
irradiance.  The  Justification  for  simply  adding  these  log- 
variances  lies  in  the  Rytov  approach  of  treating  the  log- 
irradiance  as  the  dependent  random  function.  Since  reflection' 
by  the  small  highly  reflective  target  results  in  essentially 
stairting  over  on  the  propagation  problem,  the  two  paths  are 
treated  as  independent  sequential  random  processes  and  hence 
the  variances  should  be  additive.  The  independence  of  the 
two  paths  is  reasonable,  since  the  beams,  both  being  diver- 
gent, are  not  comparable  in  size  anywhere  except  near  half- 
way between  the  transmitter  and  target.  Therefore,  the 

turbulent  effects  are  different.  This  would  be  more  rigor- 

2 

ously  correct  if  were  handled  as  a function  of  position, 
but  the  approach  should  give  reasonable  values  of  the  log- 
variance  at  the  receiver. 

Note  that  the  distributions  are  calculated  and  plotted 
for  the  irradiance  on  axis  at  the  receiver.  Beam  wander 
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effects  are  not  Included  for  the  return  path  and  beam 
receiver  interaction.  After  RINPIH  is  complete.  TARGET  is 
called  again,  using  the  average  irradiance  and  the  highest 
and  lowest  expected  power  in  the  presence  of  wander  (cal- 
culated by  TURB.  along  with  URPP  and  BRPP.  by  the  same 
method  as  UPRLIM  and  BTMLIM).  If  beam  wander  is  significant 
(JERKY  not  equal  to  zero),  a statement  is  printed  out 
including  the  average  output,  the  high  and  low  expected  out- 
put. These  high  and  low  values  correspond  to  well  out  on 
the  wings  of  the  distribution,  since  they  are  away  from  the 
average  value  by  one  variance  on  the  forward  path,  and  then 
another  (different)  variance  on  the  reflected  path.  Hence, 
the  probability  that  a value  would  fall  outside  this  range 
should  be  less  than  10%.  but  not  negligible.  If  wander  is 
negligible  (JERKY  equals  zero),  only  the  average  value  of 
the  output  is  printed,  and  a direct  ratio  between  the  axial 
izradiances  and  the  powers  can  be  assumed. 

No  units  are  specified  with  the  output  printout.  This 
is  because  the  receiver  sensitivity  is  used  in  the  place  of 
the  reflectivity.  The  units  of  the  answer  for  output  depend 
upon  the  sensitivity  value  used.  The  default  value  is  one. 
so  that  the  output  is  in  watts  of  optical  power  on  the 
receiver. 
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VI,  Combined  Model 

The  various  atmospheric  effects  are  combined  In  a pro- 
gram called  PULSE.  This  program  handles  data  Input  and 
parts  of  the  written  output.  All  of  the  calculations  of  the 
various  parameters  are  performed  In  subroutines.  The  major 
subroutines  are  described  below. 

Input 

The  Inputs  required  for  the  program  are  on  formatted 
cards.  All  of  the  Input  parameters  have  default  values  and 
input  cards  change  these  values.  The  various  cards  and 
their  required  entries  are  below.  The  formats  required  are 
listed  after  the  card  description.  The  name  of  the  card 
starts  In  colxann  one  and  Is  not  followed  by  any  punctuation. 
Notice  that  all  numerical  fields  start  In  columns  10,  20, 

30,  etc.  The  given  formats  start  In  column  10. 

NUMBER  requires  an  Integer  1,  2 or  3 in  column  10. 

FORMAT  (II) 

BEAM  sets  the  values  of  wavelength  (microns) , 

diameter  (meters),  divergence  (radians), 
peak  power  (watts),  pulse  duration  (seconds) 
and  pulse  repetition  time  (seconds).  These 
are  the  beam  parameters  for  beam  number  3. 
Beams  1 and  2 remain  set  at  default  values 
of  the  1.06  and  .53  micron  beams. 

FORMAT  (7(E  7,1,  3x)) 


r 

RANGE 

TARGET 

TURB 

SETTING 

WEATHER 

RECEIVER 

NOPLOT 

PRINT 

ENDPROB 


transmitter  to  target  (kilometers). 

FORMAT  (F  10.3) 

target  diameter  (meters),  reflectivity, 
divergence  of  the  be2un  after  reflection 
(radians).  FORMAT  (3(E  7.1,  3x)) 

effective  refractive-index  structure  con- 
stant  (CNSQ)  (m““),  transverse  wind  veloc- 
ity (m/sec).  FORMAT  (E  7.1,  3x,  F 10.3) 

requires  one  of  the  following  words  begin- 
ning in  column  10.  URBAN,  RURAL,  MARITIME 
or  NCM^STAND.  If  NONSTAND  is  entered,  curve 
fitting  parameters  C and  D are  required  in 
columns  30  and  40.  FORMAT  (AS,  llx,  2F  10.4) 

temperature  (°C) , atmospheric  pressure  (mbar) , 
visibility  (km),  dew  point  (°C). 

FORMAT  (4(E  7.1,  3x)) 

receiver  diameter  (meters),  receiver  sensi- 
tivity (appropriate  units).  FORMAT  (2(E  7.1, 
3x)) 

no  numerical  entry  required.  This  suppresses 
the  plots  of  irradiance  probability  distri- 
bution. 

no  numerical  entry  required.  This  prints 
output  of  intermediate  results  and  comments. 
Primarily  useful  if  problems  are  encountered 
in  the  program. 

no  numerical  entry  required.  This  tells  the 
program  that  no  more  input  cards  are  to  be 
read. 


1 
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The  default  values  set  in  the  program  are  equivalent  to 


an  input  deck  ofi 


NUMBER 

2 

RANGE 

8.0 

TARGET 

0.025 

1.0 

0.003 

TURB 

5.0  E-15 

2.0 

SETTING 

URBAN 

WEATHER 

25.0 

975.0 

15.0 

RECEIVER 

0.0762 

1.0 

ENDPROB 

If  no  ENDPROB  is  included,  the  end-of-file  will  perform 
the  same  function.  The  NUMBER  card  sets  the  number  of  beams 
to  be  considered.  As  currently  set  up,  only  beam  3 Cein  be 
varied.  This  is  done  in  order  to  be  able  to  set  the  values 
for  the  variable  beam,  while  using  the  fixed  values  for  the 
other  two  beams.  It  is  not  easy  to  investigate  the  results 
of  changing  beam  parameters.  A few  minor  changes  to  the 
program,  as  listed  in  Appendix  £,  allow  the  inclusion  of  a 
NUMBER  card,  followed  by  the  appropriate  ntunber  of  BEAM 
cards  (up  to  3),  to  change  all  the  beam  parameters,  rather 
than  Just  beam  3. 

Subroutines 

The  major  subroutines  are  listed  below,  with  a summary 
of  tdiat  each  one  does.  The  basis  of  each  calculation  is 
discussed  in  the  appropriate  chapter. 

AFCRL  sets  the  value  of  molecular  absorption 

coefficient  (KAPPAM)  equal  to  .001  and 
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AROSOL 


RAYLEE 


TURB 


AXIAL 


TARGET 


RTNPTH 


flags  wavelengths  where  absorption  may 
require  better  calculations. 

calculates  the  aerosol  extinction  coeffi- 
cient (SIGMA)  based  upon  approximations  to 
the  various  aerosol  models. 

calculates  the  molecular  scattering  coeffi- 
cient (SIGMAM)  based  upon  Rayleigh  scat- 
tering. The  main  program  then  sums  the 
three  coefficients  to  give  the  total 
extinction  coefficient  (SIGMA),  This  Is 
used  to  get  transmittance. 

calculates  the  basic  Irradlance  scintillation 
parameters  by  use  of  various  formulas  and 
function  subroutines.  It  also  calculates 
characteristic  scintillation  time,  SCNTIM. 

calculates  the  average  Irradlance  and  the 
beam  wander  parameters.  Including  the 
characteristic  wander  time,  WNDTIM. 

determines  the  amount  of  power  reflected  by 
the  target,  based  upon  Initial  power,  beam 
parameters,  transmittance,  target  size, 
reflectivity  and  the  turbulence  character- 
istics. Both  characteristic  times  are 
Involved  In  determining  the  pulse -to -pulse 
behavior  of  reflected  power.  It  Is  used 
again  to  determine  output  of  the  receiver. 

calculates  the  total  power  reaching  the 
receiver  plane  and  the  Irradlance  distri- 
bution. It  calls  TURB  to  evaluate  the 
effects  of  txirbulence.  It  then  calls  PROB 
or  PROBZ. 
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PROS 


] 

calculates  euid  plots  (if  required)  the 
Irradlance  probability  distribution  at  the 
receiver  based  upon  the  two -parameter  log- 
normal distribution. 

PROBZ  calculates  and  plots  (If  required)  the 

Irradlance  probability  distribution  at  the 
receiver  based  upon  the  three -parameter 
truncated  log-normal  distribution. 

Output 

The  output  Is  basically  In  two  forms,  printed  and 
plotted.  The  printed  output  (Fig.  13)  first  repeats  the 
Input  cards.  It  lists,  for  each  beam,  the  beam  parameters, 
turbulence  structure  parameter  and  transverse  wind  speed. 

It  then  tells  which  probability  distribution  Is  assumed  for 
the  given  parameters,  followed  by  the  mean,  median  and  upper 
and  lower  quartlle  points  of  the  Irradlance  distribution. 

The  next  statement  addresses  the  effect  of  beam  wander  at 
the  target.  If  the  parameter  JERKY  was  set  equal  to  zero  In 
the  wander  calculations,  the  phrase  "wander  Is  not  signifi- 
cant” Is  printed.  If  JERKY  Is  not  equal  to  zezro,  the 
statement  "Distribution  Is  modified  by  beam  wander"  Is 
printed.  This  serves  as  a warning  that  there  has  been  a 
wander  factor  added  to  the  log-variance  and  that  the  form 

i 

of  the  distribution  Is  probably  not  accurate.  With  diver- 
gent beams,  It  Is  unlikely  that  this  situation  occurs. 

The  last  two  statements  concern  effects  at  the  receiver. 

The  first  Is  a statement  derived  from  the  relationships 

between  the  ctiaracterlstlc  times  of  the  pulsed  beam,  beam  | 
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wander  and  scintillation.  The  effects  of  wander  and  scin- 
tillation are  each  described  as  random,  varying  slowly  or 
averaged  out.  If  a characteristic  time  (WNDTIM  or  SCNIIM) 
is  less  than  the  pulse  duration  (very  unlikely),  the  effect 
Is  averaged  out.  since  each  pulse  sees  some  sort  of  average 
effect.  If  a characteristic  time  is  longer  than  the  pulse 
repetition  time,  the  effect  varies  slowly,  with  successive 
pulses  seeing  a slightly  different  situation.  If  a char- 
acteristic time  is  between  the  pulse  duration  and  pulse 
xrepetition  time  (the  most  likely  situation),  the  effect  is 
random,  since  each  pulse  sees  a changed  situation.  Although 
this  is  calculated  for  the  return  path,  the  parameters  are 

f,- 

close  enough  that  both  paths  most  likely  show  the  same 
behavior.  The  last  statement,  as  described  in  Chapter  V, 
gives  the  output  from  the  receiver. 

The  plotted  output  (Fig.  14)  gives  a plot  of  the  assumed 
probability  distribution.  The  curves  are  computed  by  mul- 
tiplying the  probability  density  at  each  point  by  an  inten- 
sity interval.  The  intervals  are  equal  in  length  along  the 
logarithm's  axis  and  centered  on  each  point.  Every  second 
point  is  marked  with  a symbol.  The  header  gives  the  type 
of  distribution,  log-normal  or  modified,  where  modified 
means  the  three -parameter  distribution.  Sigma  is  the  square 
root  of  the  variance  used  to  calculate  the  curve.  The  wave- 
length and  base  10  log  of  CNSQ  are  printed  to  help  identify 
the  plot. 


70 


I 


j 

I 

I 

An  additional  pirinted  output  is  available.  This  is  a 
pzrintout  of  all  entries  in  the  prograun  written  to  XAPE7.  An 
example  is  given  in  Appendix  B for  a single  beam  problem. 

The  numerical  values  of  selected  intermediate  results  are 
printed  out.  Some  of  these  may  be  of  interest  in  themselves* 
such  as  the  scattering  and  absorption  coefficients  or  the 
I vacuum  and  wander  induced  beam  spread.  Others  are  useful 

only  for  checking  on  problems  in  the  program*  such  as  the 
number  of  terms  used  in  the  hypergeometric  series. 

Results 

It  can  be  seen  from  the  samples  included  here*  which 

should  correspond  to  a sunny  day  with  moderate  turbulence* 

that  the  probable  range  of  irradiance  at  the  receiver  varies 

over  several  orders  of  magnitude.  Furthermore*  the  temporal 

behavior  is  such  that  the  irradiance  varies  reindomly  from 

pulse -to -pulse.  Hence*  even  in  this  fixed  aim*  fixed  target 

situation*  it  is  clear  that  successive  pulses  returned  from 

the  target  may  be  expected  to  differ  by  as  much  as  two 

orders  of  magnitude.  Indeed*  from  Fig.  13*  it  can  be  seen 

for  the  1.06  micron  beam  that  only  half  of  the  pulses  fall 

within  a range  which  goes  from  3.28  x 10"^  to  1,97  x 10"^ 

2 

w/m  * a factor  of  6.0. 

Some  early  results  of  experimental  measurements  (Ref  27) 
are  plotted  in  Fig.  15.  The  units  on  the  abscissa  are  essen- 
tially power  received  per  pulse*  although  it  is  plotted  in 
terms  of  an  effective  target  cross-section.  This  rather 
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arbitrary  scale  avoids  the  problems  of  calculation  involved 
in  obtaining  exact  values,  and  is  appropriate  for  the  use  on 
a fixed  range,  such  as  is  involved  here.  The  ordinate  is  a 
linear  plot  of  number  of  pulses  received,  plotted  as  a per- 
centage of  total  pulses.  The  plot  does  not  go  to  lower 
values  of  the  received  power  because  pulses  in  that  range 
fall  into  the  noise  of  the  receiver,  so  that  an  extremely 
large  nxjmber  of  pulses  is  counted  in  a small  interval  of 
powers.  Virtually  all  of  the  left  side  of  the  distribution 
ends  up  being  counted  at  the  level  represented  by  the  noise. 

From  the  actual  tabulated  data,  the  upper  and  lower 
qxiartile  points  are  found  to  be  about  1,1  x 10^  and  3,3  x 10^, 
a factor  of  3,3,  This  corresponds  fairly  well  to  the  theo- 
retical factor  of  6,0,  especially  since  both  the  theoretical 
and  experimental  data  are  for  "typical  days,"  with  no 
attempt  to  actually  match  the  weather  and  ttirbulence  con- 
ditions. 

Several  more  examples  of  the  theozretical  output  are 


included  in  Appendix  B 


i 


' } 

■ 

VII.  Conclusions  I Discussion  and  Recommendations 

Conclusions 

The  model  as  developed  here  seems  to  fit  fairly  closely 
the  data  obtained  experimentally,  at  least  at  1.06  microns. 

The  comparison  made  in  Chapter  VI  was  a very  rough  one.  but 

I was  sufficient  to  show  a general  agreement  in  form.  It  is 

clear  that,  on  days  with  fairly  good  visibility  so  that 
aerosol  losses  are  not  significant,  the  dominant  effect  on 
the  received  power  will  be  atmosi^eric  turbulence.  On  a 
normal  sunny  day,  the  power  received  can  be  expected  to  vary 
randomly  on  a pulse -to -pulse  basis.  This  variation  can 
easily  span  two  orders  of  magnitude  of  received  power. 

I Discussion 

I A system  employing  pulsed  laser-beams  such  as  these  in 

! 

||  any  sort  of  tracking  or  aiming  problem  must  be  able  to  deal 

with  pulse -to -pulse  variations  of  two  orders  of  magnitude 

S 

j;  even  when  aiming  is  exact.  Hence,  the  tracking  algorithms 

I must  take  this  into  account.  Furthermore,  it  would  seem  j 

that  such  a system  would  be  limited  to  situations  idiere  the 

target  is  much  more  reflective  than  the  background,  since 

nearby  reflections  even  two  orders  of  magnitude  less  might 

be  mistaken  for  target  returns.  For  a good  quality  comer 

reflection,  this  should  present  no  difficulties  in  most  j 

cases.  If  the  use  of  less  ideal  reflection  is  anticipated. 
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a study  of  this  problem  must  be  Included  in  the  analysis  of 
the  tracking  problem. 

The  agreement  between  the  theory  and  experiment  in  the 
.53  micron  case  is  not  nearly  as  good  as  in  the  1.06  micron 
case  (Ref  27).  This  has  only  recently  been  noted,  and  no 
effort  has  been  made  here  to  investigate  the  problem.  A 
starting  point  would  be  to  agree  on  units  for  the  power  (or 
irradiance)  level  and  intervals,  so  that  it  would  be  possible 
to  check  the  actual  levels  of  power  received  against  the 
theoretical  prediction.  This  is  especially  important  for 
the  .53  micron  beam,  since  the  whole  distribution  is  shifted 
to  lower  power  levels,  closer  to  the  system  noise  level,  and 
it  is  therefore  important  to  know  what  part  of  the  distri- 
bution is  being  detected. 


Recommendation  s 

Further  study  needs  to  be  done  to  explain  the  .53  micron 
results,  before  the  model  can  be  used  with  any  degree  of 
I assurance  as  to  its  accuracy.  Furthermore,  there  needs  to  be 

a coordinated  effort  to  compare  weather  observations  (meteoro- 
logical measurements  as  well  as  subjective  descriptions  of  the* 

2 

air  mass  condition)  with  values  of  obtained  by  fitting  the 
curves  to  the  data.  That  has  not  been  done  here. 

If  Davidson  and  Gonzalez -del -Valle  (Ref  9).  or  others, 
improve  the  probability  distribution  description  for  inter- 
mediate and  strong  turbulence,  the  distribution  functions 
Included  here  need  to  be  improved  as  well.  In  particular, 
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no  special  provision  has  been  made  here  for  the  case  of 
Oj  greater  than  100.  For  this  caset  it  may  be  appropriate 
to  include  a Rayleigh  distribution  (Ref  10t81-82).  An 
explicit  theoretical  form  for  the  parameters  of  such  a 
distribution  has  not  been  derived • however.  Current  work 
(Ref  1)  is  proceeding  in  many  places  on  understanding  these 
probability  distributions,  and  any  systems  designer  needing 
accurate  values  of  the  expected  power  fluctuations  needs  to 
keep  abreast  of  the  current  efforts. 
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Appendix  A 
Program  Listing 


This  appendix  contains  a complete  listing  of  the  main 
program  and  subroutines.  Liberal  use  has  been  made  of 
comments  throughout.  The  function  of  each  subroutine  has 
been  described  in  the  appropriate  chapter. 
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Appendix  B 
Sample  Results 


This  appendix  contains  more  sample  output.  The  first 
type  shown  is  the  expanded  printed  output  obtained  by  the 
PRINT  card  (pp,  121-130),  Next  is  the  basic  output  for  a 
fairly  strong  turbulence  case,  = 1.0  x 10"^^,  followed  by 
a case  of  weak  turbulence,  where  = 1,0  x 10”  , Both  of 
these  include  the  pxrinted  and  the  plotted  results. 
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The  Air  Force  Geophysical  Laboratory  has  been  studying 
molecular  absorption  in  the  atmosphere  for  many  years.  The 
basic  approach  has  been  to  determine  individual  absorption 
lines  for  each  molecule  and  each  set  of  energy  levels.  This 
has  been  accomplished  both  by  theoretical  calculations  amd 
by  laboratory  measurements.  Although  only  the  seven  major 
constituents  of  the  atmosphere  (in  terms  of  absorption) 
have  been  considered,  over  150.000  individual  absorption 
lines  have  been  tabulated  (Ref  21).  This  data  base  is 
constantly  updated  as  new  or  improved  information  becomes 
available.  It  is  contemplated  to  either  add  to  this  com- 
pilation. or  compile  separately,  a similar  collection  for 
the  major  pollutants,  notably  hydrocarbons  (Ref  3),  This 
would  be  of  particular  interest  in  the  visible  and  near-lR. 
since  in  most  of  this  range  little  absorption  is  present 
due  to  the  primary  constituents. 

The  compilation  is  on  magnetic  tape  and  consists,  in 
addition  to  identification  by  molecule,  isotope  and  quantum 
levels  of  the  transition  involved,  the  following  information i 

V ■ frequency  (wavenumber)  in  cm*^ 

S ■ line  intensity  in  cm’Vnolecule  cm"^  at  296®K 
a ■ half -width  in  cm”^atm  at  296 °K 

energy  of  the  lower  state  in  cm~^ 


A great  deal  of  effort  has  gone  into  understanding  the 
line- shape  of  these  transitions  and  into  modeling  the  effects 
of  temperature  and  pressure  on  the  line shapes. 

The  data  base  and  detailed  knowledge  of  the  lineshapes 
have  led  to  the  development  of  four  different  types  of  com- 

r ' «. 

puter  cards  to  extract  useful  information  from  the  data.  An 
example  of  each  is  described  below. 

Program  LIN 

Program  LIN  is  the  basic  brute  force  approach  to 
obtaining  high  resolution  absorption  spectra.  It  requires 
that  the  number  densities  of  the  seven  absorbing  species  be 
input  at  up  to  15  pressure  levels  of  interest,  as  well  as 
the  pressure  and  temperature  at  each  level.  Then  the  limits 
of  the  frequency  range  to  be  considered,  Vl  and  V2,  eind 
the  interval  at  which  results  are  printed  (or  more  commonly, 
plotted),  DELV,  are  entered.  Additional  parameters  needed 
are  DV,  the  interval  at  which  monochromatic  calculations 
are  madei  BOUND,  the  frequency  range  from  the  frequency  of 
calculation  to  which  the  program  searches  for  lines  and  A, 
the  half -width  of  a triangular  slit  function  which  is 
convolved  with  the  monochromatic  results.  Another  parameter, 
DEPTH,  can  be  used  to  reject  lines  which  have  an  optical 
depth  (line  strength  x molecular  number  density)  below  a 
certain  level.  The  general  slant  path  case  also  requires 
a path  description,  generally  a range  and  the  secant  of  the 
zenith  angle. 


A general  flow  chart  of  the  program  is  given  in 
Figs.  C-1  and  C-2.  No  details  of  the  calculations  are  given, 
and  large  parts  of  the  program  are  not  included  at  all. 

Those  omitted  are  essentially  bookkeeping  sections  which 
are  very  necessary  for  the  proper  operation  of  the  program. 

At  the  initial  stage  of  reading  in  the  line  parameters,  if 
the  frequency  range  from  Vl -A -BOUND  to  V2  + A + BOUND  contains 
over  3000  lines,  the  calculation  cannot  proceed  directly, 
but  must  be  broken  up  into  segments.  The  calculations  of 
tremsmittance  are  completed  for  one  segment,  up  to  a fre- 
quency W,  which  is  (A+BOUND)  below  the  upper  line  included. 
Then,  all  lines  above  W-(A+BOUND)  must  be  iretained  while 
new  lines  are  read  in  so  that  information  in  the  boundary 
areas  is  not  lost.  If  A and  BOUND  are  not  carefully  chosen, 
this  can  lead  to  adding  only  a few  lines  each  time  or  even 
failure  of  the  whole  approach.  This  is  especially  true  in 
spectral  regions  with  many  lines.  This  effect  can  be  mini- 
mized to  some  extent  by  using  higher  values  of  DEPTH,  thus 
rejecting  more  lines.  In  dense  regions  this  is  not  satis- 
factory either,  since  the  cumulative  effort  of  several  weak 
lines  may  be  significant. 

The  discussion  above  gives  a very  sketchy  idea  of  the 
line-by-line  approach  and  a small  feeling  for  the  amount  of 
calculation  involved.  Computer  processing  time  can  easily 
run  into  the  hundreds  of  seconds  for  spectral  regions  of 
reasonable  size. 
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Fig.  C-2.  Homogeneous  Path  Calculation 
(Program  LIN) 
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An  example  of  the  output  from  LIN  is  given  in  Fig.  C-3. 

i 

It  is,  of  course,  also  possible  to  get  a numerical  listing 

f 

of  the  output. 

There  are  versions  of  this  program  which  will  accept 
convolution  fxinctions  other  than  the  triangular  slit 
function.  This  would  be  useful  for  a particular  laser  line- 
shape  or  for  a receiver  spectral  response,  vrtiere  very 
accurate  results  are  required. 

Program  LASER 

Program  LASER  (Ref  24)  is  designed  to  give  strictly 
monochromatic  results.  The  calculation  of  the  transmitt£uice 
is  exactly  the  same  as  that  for  a single  frequency  in 
Program  LIN.  The  calculation  is  repeated  for  six  different 
standard  atmosphere  models,  hence  giving  a feel  for  the 
range  of  absorption  coefficient  values  to  be  expected  under 
various  conditions.  A sample  of  the  output  is  given  in 
Fig.  C-4.  (This  is  from  an  older  version  of  the  program, 
which  included  only  five  standard  models.  It  was  selected 
to  show  the  1.06  micron  results.) 

Note  that  the  output  gives  values  for  the  Rayleigh 
scattering  coefficient,  as  well  as  aerosol  attenuation 
coefficients  for  two  common  aerosol  models.  These  will  not 
be  discussed  here,  but  it  is  clear  that  this  form  of  out- 
put is  excellent  for  giving  an  estimate  of  the  relative 
importance  of  various  effects. 
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Fig.  C-4.  Results  of  Program  LASER 


Since  there  is  no  information  available  for  other  fre- 
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quencies,  the  choice  of  V»  A*  BOUND  and  DEPTH  are  critically 
important  in  this  case.  The  problem  mentioned  in  Chapter  111 
concerning  precisely  defined  frequencies  (or  wavelengths) » 
and  carefully  picking  A to  describe  the  line>width  of  the 
laser  beam  could  make  a difference  of  orders  of  magnitude  in 
the  absorption,  particularly  for  very  narrow  laser  lines. 

Program  HIRACC 

Program  HIRACC  (:Iigh  Resolution  Absorption  Coefficient 
Code)  is  a program  developed  at  AFGL  (Ref  7)  to  replace  the 
line-by-line  calculation.  The  version  described  is  the 
basic  version  for  the  uniform  path  case,  using  a Lorentzian 
line  shape,  suitable  for  low  altitude  use.  A program 
basically  the  same,  called  HIRACV,  which  uses  an  approximation 
to  the  Voight  lineshape  and  is  valid  at  all  altitudes,  is 
being  incorporated  into  a much  more  general  code  called  (at 
least  at  present)  FASCODE,  for  Fast  Atmospheric  Signatvire 
Code  (Ref  8).  This  code  will  give  a capability  to  calculate 
not  only  absorption,  but  also  atmospheric  emission,  for  the 
general  slant  path  case,  either  looking  up  or  down  through 
the  atmosphere.  This  will  give  an  essentially  complete 
characterization  of  an  atmospheric  path,  for  molecular 
effects. 

Fig.  C-5  is  taken  from  a preliminary  report  on  program 
HIRACC,  and  compares  the  results  to  those  of  a conventional 
line-by-line  calculation.  Note  particularly  the  saving  in 
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Fig.  C*5.  Coaparison  of  Output  Froa  HIRACC  (New)  and  LIN  (Conventional) 


calculation  time.  The  savings  in  time  have  been  made 
primarily  by  approximating  the  Lorentz  line  shape  by  the  sum 
of  three  functions.  The  Lorentzian  falls  to  zero  very 
slowly  in  the  far  wings,  and  this  makes  use  of  uniform 
intervals  in  the  computer  calculation  very  uneconomical  of 
time.  The  three  functions  into  which  the  Lorentz  shape  is 
decomposed  behave  as  fast,  slow  and  very  slow  functions  of 
frequency,  and  are  sampled  at  different  intervals.  This 
results  in  about  a factor  of  five  improvement  in  calculating 
time,  when  the  appropriate  sampling  intervals  are  chosen. 
Another  factor  of  two  comes  from  using  tabulated  values  for 
the  three  functions,  rather  than  calculated,  and  from  more 
optimized  selection  of  calculational  steps  in  loops  which 
are  used  many  times. 

A major  difference  between  the  approaches  of  this 
program  and  program  LIN  is  that  HIRACC  calculates  all  the 
absorption  coefficients  due  to  one  line,  then  moves  to  the 
next  line,  etc..  Program  LIN,  on  the  other  hand,  calculates 
the  contributions  from  all  lines  at  a given  frequency,  then 
indexes  in  frequency  and  repeats.  This  requires  the  HIRACC 
program  to  sum  the  absorption  at  each  frequency  to  obtain 
the  output,  but  it  reduces  the  number  of  references  to  the 
line  data.  BOUND,  in  the  case  of  HIRAAC  is  automatically 
set  to  48  times  the  half -width  of  the  line  (other  multiples 
are  possible,  if  the  decomposed  functions  are  changed). 

Lines  are  not  rejected  by  a DEPTH  parameter,  and  calculations 
of  absorption  coefficient  are  made  between  VI  and  V2,  at 
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intervals  of  DV,  with  no  slit  function  convolution.  As  In 
LIN,  computations  are  broken  into  segments  (called  Panels) 
when  necessary. 

Output  from  the  basic  program,  as  in  the  case  of  LIN, 
is  a binary  file  containing  a header  record  with  all  the 
input  parameters,  followed  by  Panel  header  information  eind 
then  absorption  coefficients  for  each  frequency  of  the  il)  ^ 

panel.  This  output  is  normally  plotted,  as  in  Fig.  4.  This 
is  a plot  of  transmittance  over  a 100  km  range  at  sea  level, 
using  the  number  densities  of  the  midlatitude  winter  model. 
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Appendix  D 

The  Log-Normai^  Distribution 


The  description  here  is  taken  from  a book  by  Aitchison 
(Ref  2)  which  primarily  addresses  the  social  sciences* 

Only  a few  of  the  mathematical  properties  of  the  distri- 
bution (or  properly*  the  distributions)  are  considered  here. 
The  basic  nature  of  the  distributions  is  that  the  logarithm 
of  the  variable,  or  the  logarithm  of  some  simple  fxinction 
of  the  variable,  is  normally  distributed. 


Two-Parameter  Log-Normal  Distribution 

The  simplest  such  distribution  is  the  two-parameter 
case.  Here,  the  two  parameters  are  (i,  the  mean  of  the  log 
of  the  variable  and  o , the  variance  in  the  log.  In  this 
case,  the  probability  density  can  be  expressed  as 


dx 


— expf. 
xo(2it)^  L 2(J^ 


(1) 


for  X > 0. 

Unlike  the  normal  distribution,  this  is  limited  to 
positive  values  of  x,  and  the  mean,  median  and  mode  axre  all 
distinct.  Fig.  D-1  is  from  Aitchison,  and  shows  a normal 
and  a log-normal  curve.  The  three  central  measures  fall 
at  the  following  values  of  x. 


151 


(Ref  2) 


8.  Normal  Distribution 


52 


Mean  » EXP(p  + %a^) 

(2) 

Mode  = EXP(p  - ah 

(3) 

Median  » £XP(p) 

(4) 

The  quantile  points  of  the  distribution 

can  most  easily 

be  found  from  those  of  the  normal  distribution,  which  are 
widely  tabulated.  If  q is  the  quantile  of  the  normal  distri- 
bution and  Q is  the  quantile  of  the  log-normal,  these  are 
related  by 

Q = EXP(n  + q<y)  (5) 

Hence,  the  lower  and  upper  quartile  points  of  the  log- 
normal are  at  x = EXP(p  - ,67a)  and  x = EXP(|i  + .67a), 
respectively.  These  are  the  values  used  in  the  program  to 
describe  the  range  in  the  printed  output. 

Three -Parameter  Log-Normal  Distribution 

The  three -parameter  function  is  basically  Just  the  two- 
parameter  distribution  with  the  variable  x displaced  by  an 
amount  T,  so  that  log(x  -T)  is  normally  distributed.  Since 
the  two-parameter  distribution  is  only  defined  for  x > 0,  the 
three -parameter  distribution  is  only  valid  for  x > T,  The 
probability  of  x <T  is  zero. 

All  of  the  central  measures  and  the  quantiles  are  Just 
displaced  by  an  amount  T,  so  that 

Mean  - T + EXP(n  + %a^)  (6) 

Mode  - EXP(4  - a^)  (7) 
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(8) 

(9) 


Median  » T + EXPCji) 

Q ■ T + EXP(|i  + qcf) 

The  reason  for  using  the  three -parameter  distribution 
(other  than  the  experimental  results  of  Ref  9)  is  that  a 
small  negative  value  of  irradiance  can  be  used  for  T, 
thus  permitting  an  irradiance  of  zero  to  be  observed • which 
should  be  admitted  as  a possibility  in  strong  turbulence. 

This  assumption  of  a negative  value  for  T,  necessarily  means 
that  a truncated  distribution  must  be  considered,  since  a 
negative  irradiance  is  physically  meaningless. 

The  truncated  distribution  is  much  more  complicated  to 
deal  with,  but  if  T is  chosen  as  a very  small  negative 
quantity,  there  should  be  little  loss  in  accuracy  in  assuming 
the  same  form  as  above  for  the  central  measures  and  quantiles, 
except  for  quantiles  which  lie  very  near  zero.  This  is 
because  only  the  effect  of  the  few  points  between  T and 
zero  are  being  neglected. 


Appendix  E 

Program  Changes  to  Vary  Beam  Parameters 


In  oztler  to  more  readily  examine  the  results  of 
changing  beam  parameters,  the  following  changes  can  be  made 
to  Program  Pulse.  The  cards  below  replace  the  same  numbered 
cards  in  the  program,  or  are  Inserted  where  Indicated  by 
numerical  order. 


NCK  = 0 000721 

2 IF(NCK.EQ.O)  CK)  TO  23  000880 

NCK  = NCK  + 1 000890 

DECODE  (70, 103,KRUD)WVLN(NCK) ,DIVERG(NCK) , 

POWER(NCK) , 000900 

-C  PD(NCK),PRT(NCK)  000910 

GO  TO  21  000912 

23  PRINT(6,24)NUMBER  000914 

24  FORMATCO  more  LINES  ^VERE  INPUT  THAN  YOU 

CALLED  FOR."  000916 

C "ALL  LINES  PAST  ", I 2, "WERE  IGNORED.")  000920 


If  these  changes  are  made,  NUMBER  must  be  Input  before 
BEAM  cards.  If  NUMBER  Is  Input  as  3 and  less  than  3 BEAM 
cards  are  Input,  the  program  will  not  run  to  completion, 
since  no  default  values  are  Included  for  beam  3.  An  addi- 
tional DATA  statement  with  beam  3 default  parameters  would 
allow  this  sort  of  Input. 
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